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INTRODUCTION 


1.1  BACKGROUND 

Omega  Navigation  System  signals  have  the  advantage  of  long  range  when  compared  to 
other  terrestrial-based  radionavigation  systems.  Because  of  their  long  range,  Omega  signals 
propagating  between  the  earth  and  the  ionosphere  encounter  many  different  perturbing  effects. 
These  effects  are  due  to  ground  conductivity  (which  vanes  over  the  earth  by  over  five  orders  of 
magnitude),  the  earth’s  magnetic  field,  conductivity  of  the  ionosphere  which  changes  with  time, 
and  direction  of  signal  propagation.  Because  of  the  associated  variations  in  signal  phase  caused 
by  these  effects,  acceptable  Omega  navigation  accuracy  is  achieved  only  if  the  signal  is 
corrected,  i.e.,  use  of  a  priori  knowledge  and  compensation  for  these  effects.  Such  a  priori 
knowledge,  or  prediction,  is  possible  (though  not  simple)  if  the  signal,  as  a  function  of  distance 
and  time,  is  sufficiently  well-behaved.  In  the  context  of  the  most  commonly-used  signal 
propagation  model  (the  waveguide-mode  model),  “sufficiently  well-behaved”  means  the  signal 
is  primarily  composed  of  a  single  “mode”  (specifically,  Mode  1  ;  see  Chapter  3  for  further 
discussion  of  signal  modes).  For  those  regions/times  in  which  the  propagated  Omega  signal 
comprises  multiple  modes  (of  roughly  equal  strength)  or  is  dominated  by  a  single  higher-order 
mode,  usable  prediction  of  the  signal  phase  is  simply  not  possible,  due  both  to  the  theoretical 
complexity  and  the  inherent  signal  instability  under  such  conditions.  Stated  another  way,  the 
phase  of  an  Omega  signal  in  the  “far  field"  (i.e.,  greater  than  about  1000  km  from  the 
transmitting  source)  varies  approximately  in  a  linear  fashion  with  distance  from  the 
transmitting  source.  Small  departures  from  this  linear  behavior  (i.e.,  Mode  1  conditions)  can  be 
predicted  but  large,  rapid  variations  and  sudden  shifts  (due  to  the  presence  of  higher-order 
modes)  cannot  be  predicted. 

Omega  propagation  corrections  (PPCs)  are  predicted  phase  values  which,  when  applied 
to  the  received  Omega  signal  phase  measurement,  provide  a  “nominal”  linear  phase  versus 
distance  relationship.  The  nominal  “value”  is  simply  the  ratio  of  the  cumulative  “idealized” 
phase  developed  by  a  signal  to  the  distance  over  which  the  signal  is  propagated  (in  general,  this 
ratio  is  called  the  wave  number).  The  numerical  value  of  the  nominal  is  not  critical;  it  is  only  im¬ 
portant  that  the  value  be  “near”  the  average  over  all  time  and  space  conditions.  The  central 
problem  in  developing  Omega  PPCs,  then,  is  to  calculate  the  variation  of  Omega  signal  phase 
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relative  to  a  reference,  i.e.,  the  nominal  phase.  Since  a  number  of  geophysical  parameters  are 
involved,  the  model  for  phase  variation  is  conveniently  partitioned  into  sub-models  which  are  ei¬ 
ther  time-dependent  (temporal)  or  space-dependent  (spatial).  The  principal  temporal 
sub-model  is  the  diurnal  sub-model  which  mainly  treats  the  solar  zenith  angle  dependence  of  the 
phase  variation.  The  spatial  sub-models  describe  effects  due  to:  geomagnetic  bearing/latitude, 
ground  conductivity,  auroral  zone,  and  polar  region.  Each  sub-model  includes  functional  forms 
and  associated  calibration  coefficients  for  computing  the  contribution  to  the  phase  variation. 

The  importance  of  Omega  PPCs  was  recognized  very  early  in  the  System’s  history.  J.A. 
Pierce  was  the  first  to  propose  simple  corrections  based  on  day/night  considerations  (Ref  1). 
The  Naval  Electronics  Laboratory  (now  known  as  Naval  Ocean  Systems  Center  (NOSC))  began 
an  extensive  theoretical  and  experimental  investigation  of  Omega  signal  propagation  prediction 
in  the  1960's  which  laid  the  groundwork  for  a  semi-empirical  Omega  phase  prediction  model. 
E.R.  Swanson  of  NOSC  led  the  development  effort  which  resulted  in  the  current  Omega  PPC 
model  structure  (Ref.  2).  A  model  based  cn  path  segmentation  and  waveguide-mode  theory  was 
available  in  the  late  1960’s  but  the  first  widely  disseminated  documentation  of  the  model 
appeared  in  1974  (Ref.  3)  when  the  U.f.  Coast  Guard  took  over  responsibility  for  Omega 
engineering,  operations,  and  user  support.  In  1980,  a  major  re-calibration  of  the  1974  PPC 
model  structure  was  undertaken  using  three-frequency  data  acquired  from  a  global  distribution 
of  monitor  sites  (Ref.  4).  The  only  change  to  the  model  structure  was  the  addition  of  a  long-term 
time  contr.bution  to  the  phase  variation.  This  addition  was  based  on  the  assumption  that  the 
ionosphere  is  “hardening”  (increased  ionization)  by  a  fixed  percentage  increase  each  year,  thus 
leading  to  a  long-term  exponential  increase  in  phase  variation.  The  validity  of  this  assumption  is 
the  subject  of  a  separate  investigation  and  is  not  addressed  in  this  report 

The  PPCs  based  on  the  1980  PPC  model  calibration  showed  a  significant  improvement 
over  those  of  the  1974  model.  However,  subsequent  usage  and  experimentation  revealed  some 
shortcomings  in  addition  to  those  which  may  be  due  to  the  predicted  long-term  variation.  Most 
of  the  reported  large  prediction  errors  (Ref.  5)  occur  when  the  signal  path  is  in  transition,  i.e.  a 
sunrise/sunset  terminator  lies  between  the  transmitter  and  receiver.  On  that  portion  of  the  path 
near  the  terminator,  the  current  diurnal  sub-model  predicts  phase  variations  from  anomalous 
injections  (“dumps”)  of  electrons  into  the  ionosphere  and  a  “dynamic"  contribution  based  on 
an  incorrect  characterization  of  D-region  dynamics.  In  addition,  transition  errors  may  be  due  to 
the  absence  of  transition  data  in  the  model’s  calibration  database.  Large  phase  prediction  errors 
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have  also  been  detected  on  day  and  night  transequatorial  paths  (Ref.  5).  These  large  errors  may 
indicate  the  need  for  a  separate  sub-model  invoked  in  the  equatorial  region  or  could,  perhaps, 
be  explained  as  an  improperly  modeled  feature  within  the  geomagnetic  bearing/latitude 
sub-model.  Similarly,  large  observed  errors  over  the  higher  latitude  regions  (Ref.  6)  could 
indicate  the  need  for  revision  of  the  polar/auroral  sub-models  or  may  be  due  to  the  unmodeled 
coupling  between  ground  conductivity  and  geomagnetic  field  direction  (with  respect  to  the 
direction  of  signal  propagation).  These  observed  errors  and  possible  model  deficiencies  point  to 
the  need  for  a  revision  of  the  model  structure. 

1.2  OBJECTTVE/APPROACH 

The  objective  of  the  effort  addressed  in  this  report  is  to  formulate  a  new  semi-empirical  PPC 
model  structure  which  retains  the  advantages  of  the  current  model  structure  and  remeaies  the  appar¬ 
ent  deficiencies  identified  in  section  1.1.  The  new  model  structure  should  conform  to  current  geo¬ 
physical/wave  propagation  models  as  much  as  possible  yet  not  be  overburdened  with  complex, 
detailed  calculations.  Being  semi-empirical  in  form,  the  model  should  be  data-driven  as  much  as 
possible,  especially  in  the  calculation  of  those  parameters  to  which  the  phase  prediction  accura¬ 
cy  is  especially  sensitive.  Therefore,  the  functional  forms  comprising  each  sub-model  should  be 
as  few  as  possible  (yet  retaining  essential  characteristic  behavior),  simple  in  form,  and  indepen¬ 
dent/uncorrelated.  The  sub-models  should  be  appropriately  combined  or  coupled,  and  be  con¬ 
sistent  with  physical  models  and  observed  data.  Finally,  the  model  should  be  structured  to 
permit  calibration  with  phase  data  recorded  on  paths  at  all  times. 

The  new  model  structure  should  include  sub-models  for  characterizing: 

1)  Time-dependent  ionospheric  parameters 

2)  Vector  geomagnetic  field 

3)  Direction  of  signal  propagation 

and  other  geophysical  quantities,  if  significantly  coupled  to  the  three  parameters  above. 
Revision  of  the  current  sub-models  characterizing  the  phase  variation  in  the  pclar/auroral  zones 
as  well  as  long-term  phase  variation  are  not  addressed  in  this  report. 

1.3  Report  Overview 

Presentation  of  the  new  model  structure  is  organized  to  first  address  the  development  of 
each  sub-model  separately  and  then  address  the  integration  of  the  sub-models  into  an  c.erall 
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structure.  This  presentation  structure  is  employed  because  the  composition  of  each  sub-model 
involves  specific/unique  concepts. 

Chapter  2  presents  the  diurnal  sub-model  which  includes  models  of  the  daytime  and 
nighttime  phase  variation.  Parameters  defining  the  important  day/night  time  “boundary”  are 
presented  and  the  predicted  boundary  times  are  compared  with  observed  data.  Since  this  chap¬ 
ter  presents  phase  variation  forms  which  are  substantially  different  from  those  of  the  current 
(1980)  model,  the  development  of  these  new  functional  forms  is  presented  in  some  depth.  Con¬ 
sequently,  in  order  to  provide  the  basis  for  understanding  the  origins  and  rationale  for  the  new 
functional  forms,  Chapter  2  includes  a  brief  review  of  the  relevant  principles  of  ionospheric 
physics 

The  single  spatial  sub-model,  which  combines  geomagnetic  field,  wave  propagation 
direction,  and  ground  conductivity  effects  is  presented  in  Chapter  3.  The  integration  of  these 
geophysical  quantities  into  one  sub-model  reflects  the  parameter  coupling  implied  by  standard 
wave  propagation  models.  The  chapter  also  presents  alternative  descriptions  of  phase  behavior 
in  a  certain  range  of  wave  propagation  directions  near  the  geomagnetic  equator. 

Chapter  4  provides  a  description  how  the  two  sub-models  are  combined  to  provide  the 
total  phase  variation/PPC.  The  temporal  forms  are  recast  as  “evolution  functions”  which  dictate 
the  changes  in  the  spatial  functional  forms  with  changing  solar  zenith  angle.  Model  calibration 
is  a  three-step  cycle:  determine  the  calibration  coefficients  for  daytime,  nighttime,  and 
parameters  for  the  day/night  time  boundary  (onset/recovery  times).  A  method  for  calibrating 
the  onset/recovery  time  parameters  is  also  presented.  Chapter  4  concludes  with  a  comparison  of 
the  current  (1980)  and  the  new  (1990)  PPC  models.  The  last  chapter  includes  a  summary,  to¬ 
gether  with  conclusions  and  recommendations. 

The  appendices  include  mathematical  details  of  some  of  the  quantitative  featu¬ 
res/processes  of  ionospheric  dynamics  needed  for  prediction  of  the  phase  variation.  Also 
included  as  appendix  material  is  a  description  of  the  analog  data  used  to  develop/test  the  new 
diurnal  sub-model. 
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2. 


DIURNAL  SUB  MODEL 


2.1  INTRODUCTION 

This  chapter  presents  tne  temporal  characteristics  of  the  predicted  Omega  signal  phase. 
These  temporal  characteristics  are  almost  entirely  manifested  by  changes  in  the  ionosphere 
(ground  conductivity  and  other  electromagnetic  properties  affecting  wave  propagation  vary 
little  with  time  at  Omega  frequencies).  Thus,  a  description  of  signal  phase  variation  over  time 
scales  of  interest  to  the  Omega  user  (more  than  a  few  minutes  but  less  than  a  year)  must  focus 
on  the  dynamics  of  the  ionosphere  in  addition  to  the  characteristics  of  wave  propagation. 

Prediction  of  signal  phase  behavior  over  long  (>650  km)  sub-ionospheric  paths  requires 
models  of  wave/ionosphere  interaction,  wave/earth-surface  interaction,  and  wave  propagation 
between  the  earth’s  surface  and  ionosphere.  As  noted  above,  the  temporal  component  of  signal 
phase  prediction  arises  from  the  changing  ionosphere  and  the  resulting  change  in  wave 
propagation  characteristics.  Wave  propagation  characteristics  are  assumed  to  follow  from  a 
waveguide-mode  model  in  which  the  path  between  transmitter  and  receiver  is  partitioned  into 
segments,  each  of  which  are  sufficiently  short  so  that  earth’s  surface  and  ionosphere  properties 
are  constant  within  the  segment.  Each  segment  can  be  viewed  as  a  section  of  an  indefinitely  long 
waveguide  having  the  same  boundary  properties  as  the  segment.  The  long  waveguide  has  a 
characteristic  phase  versus  distance  relationship  which  differs  from  that  of  free  space.  This 
difference,  referenced  to  a  path  segment  length,  is  the  incremental  phase  contributed  by  the 
path  segment.  The  sum  (over  all  path  segments)  of  the  incremental  phase  contributions  is  added 
to  the  total  free-space  phase  change  along  the  path  to  obtain  the  total  predicted  path  phase. 

At  the  local  path  segment  level,  the  ionosphere  is  predominantly  controlled  by  direct 
solar  radiation  from  the  point  at  which  the  sun  is  directly  overhead  to  a  threshold  solar  zenith 
angle  (see  Fig.  2.1-l(a)).  This  range  of  solar  zenith  angle,  when  referenced  to  a  path  segment, 
refers  to  an  illumination  condition  known  as  local  day.  For  solar  zenith  angles  larger  than  the 
threshold  value,  the  path  segment  ionosphere  is  controlled  by  ionization  sources  other  than 
direct  solar  radiation  and  the  path  segment  illumination  condition  is  defined  as  local  night. 

A  path  is  composed  of  numerous  segments  and,  as  a  result,  its  illumination  conditions 
differ  from  those  of  a  single  segment.  Path  night  is  said  to  occur  when  all  path  segments  are  in  a 
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b)  Path  Temporal  Conditions 

Figure  2.1-1  Local  (Path  Segment)  and  Path  Temporal  Conditions 


local  night  condition.  The  path  darkness  fraction  (PDF)  is  the  number  of  path  segments  in  a  local 
night  condition  divided  by  the  total  number  of  segments  (in  which  the  wave  interacts  with  the 
ionosphere).  Thus,  for  path  night,  PDF=1.0.  Path  sunrise  is  defined  as  the  condition  in  which  the 
PDF  decreases  in  time  (i.e.,  path  segments  are  changing  from  a  local  night  to  a  local  day 
condition).  Path  day  occurs  when  PDF=0.0,  i.e.,  when  all  path  segments  are  in  a  local  day 
condition.  Path  sunset  is  said  to  occur  when  the  PDF  increases  on  the  path.  Figure  2.1-1 
illustrates  and  compares  local  and  path  temporal  conditions. 


Once  phase  behavior  at  the  local  day  and  local  night  ionosphere  has  been  established,  it 
is  critical  to  determine  the  time  at  which  each  of  these  conditions  begins  and  ends.  The  first 
arrow  shown  (following  T=0)  in  Fig.  2.  l-l(b)  indicates  the  time  at  which  the  threshold  solar 
zenith  angle  is  attained  at  the  easternmost  path  segment  and  is  known  as  the  sunrise  onset  time. 
The  second  arrow  indicates  the  time  at  which  the  threshold  solar  zenith  angle  is  attained  at  the 


westernmost  path  segment  and  is  known  as  the  sunrise  recovery  time.  Similarly,  the  third  arrow 
indicates  the  sunset  onset  time  and  the  fourth  arrow  the  sunset  recovery  time. 


The  importance  of  accurately  specifying  onset/recovery  times  is  illustrated  in  Fig.  2.1-2. 

From  the  figure,  it  is  clear  that  the  phase  error  (difference  between  observed  and  predicted 

phase)  during  path  sunrise  is  a  monotonically  increasing  function  of  both  the  onset-time 

prediction  error  and  the  slope  of  the  sunrise  phase  decrease.  The  relationship  between  phase 

error  and  onset  time  error  is  given  approximately  by 

A  ^  _  0.03LATo 
- — - 

A<I> 

where  L  is  the  path  length  (in  Mm),  AT0  is  the  error  in  onset-time  prediction  (in  minutes),  Ad>  is 
the  longitude  difference  between  the  path  end-points  (in  degrees),  and  A0e  is  the  resulting 
phase  prediction  error  (in  cycles).  For  a  given  onset-time  prediction  error,  the  phase  error  is 
seen  to  be  larger  for  longer  paths  and  those  that  are  more  north-south  oriented  (A<l>  smaller). 


PREDICTION  ERROR 


Figure  2.1-2  Phase  Prediction  Error  due  to  Onset  Time  Prediction  Error 


2-3 


This  chapter  is  organized  to  address  the  three  most  important  aspects  of  diurnal  phase 
predication  noted  above: 

•  Signal  phase  perturbations  due  to  a  local  daytime  ionosphere 

•  Signal  phase  perturbations  due  to  a  local  nighttime  ionosphere 

•  Parameters  describing  path  sunrise/sunset  onset/recovery  times 

Since  the  description  of  phase  behavior  under  daytime/nighttime  ionosphere  requires  some 
familiarity  with  the  dynamics  and  morphology  of  the  lower  ionosphere,  the  chapter  begins  with 
a  brief  background  in  ionospheric  physics. 


2.2  BACKGROUND:  IONOSPHERIC  PHYSICS 
2.2.1  Ionospheric  Dynamics 

In  simplest  terms,  the  ionosphere  is  that  portion  of  the  upper  atmosphere  (above  about 
70  km)  in  which  a  small  but  important  component  of  the  particle  population  consists  of  charged 
particles  whose  numbers  are  maintained  as  a  result  of  the  competing  forces  of  solar 
photoionization  and  reconstitution  into  neutral  atoms/molecules  as  a  result  of  various 
mechanisms.  For  a  given  particle  species,  this  balance  is  described  by  a  continuity  equation 
which  states  that  (for  a  given  volume): 

Rate  of  change  of  charged  particle  number  =  (charged  particle  production  rate)  — 

(charged  particle  loss  rate) 

The  types  of  charged  particles  found  in  the  ionosphere  include  electrons,  negative  ions  and 
positive  ions.  The  specific  components  which  are  ionized  in  the  lower  ionosphere  are  mainly 
atomic/molecular  forms  of  nitrogen,  oxygen,  and  nitric  oxide. 

In  tracing  the  balance  of  positive  ions  (other  species  give  similar  results),  it  is  seen  that 
positive  ions  are  lost  when  combined  with  electrons  or  negative  ions.  The  rate  of  loss  is 
proportional  to  both  the  positive  ion  density  and  the  combining  species  density.  Thus,  the 
continuity  equation  becomes  (Ref.  7) 

=  q  -  OcN  +  Nc  -  cijN + N_  (2.2-1) 

dt 

where  N+  is  the  positive  ion  density,  Nc  is  the  electron  density,  N_  is  the  negative  ion  density,  q 
is  the  production  function,  ae  is  the  electron-positive  ion  recombination  coefficient,  a\  is  the 
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positive  ion-negative  ion  recombination  coefficient,  and  t  is  the  time.  Since  charge  neutrality 
holds  for  time  scales  of  minutes  throughout  the  ionosphere, 

N+  =  N_  +  Ne  =  (1  +  A)NC  (2.2-2) 

where  A  =  N_/Ne  is  the  negative  ion  ratio.  Since  A  is  practically  independent  of  time, 
substitution  of  Eq.  2.2-2  into  Eq.  2.2-1  yields 

=  q'  -  akNi  (2.2-3) 

dt 

where  q'  =  q/(l  +  A)  and  cle'  ~  at  +  ka\  is  known  as  the  effective  recombination  coefficient. 
Equation  2.2-3  can  be  solved  analytically  for  very  simple  forms  of  q';  otherwise  approximation 
methods  or  numerical  integration  must  be  employed. 

2.2.2  Wave  Reflection  Height 

The  lower  ionosphere  is  frequently  defined  as  a  “weakly  ionized  plasma”  which  means 
that  it  does  not  consist  entirely  of  charged  particles.  In  fact  only  a  small  (but  important)  fraction 
of  the  particles  in  the  lower  ionosphere  is  charged.  If  the  charged  particles  of  the  ionospheric 
plasma  are  separated  and  then  released,  the  plasma  oscillates  at  a  frequency  known  as  the 
plasma  frequency.  No  electromagnetic  wave  with  a  frequency  less  than  the  plasma  frequency  can 
propagate  through  the  plasma  (without  losing  energy  to  the  particles)  since  longer  wave  periods 
impart  larger  spatial  amplitudes  to  the  particles  which  drives  them  out  of  their  equilibrium 
states.  For  a  plasma  of  electrons  and  positive  ions,  the  plasma  frequency  (MKS  units)  is  given 
by 

(op  = 

where  Ne  is  the  electron  density,  e  is  the  electron  change,  m  is  the  mass  of  the  electron  and  e0  is 
the  permittivity  of  free  space. 

Thus,  an  electromagnetic  wave  of  frequency  to  will  propagate  through  regions  of  the 
ionosphere  for  which 

Ne 

The  electron  density  generally  increases  with  height  since  fewer  electron-ion  collisions  occur  at 
higher  altitudes  (where  total  particle  density  is  smaller)  and  thus  electron  loss  due  to 
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recombination  is  less  frequent.  Hence,  an  electromagnetic  wave  of  frequency  (o  generated  at  the 
ground  travels  upward  until  it  encounters  an  ionospheric  layer  with  electron  density 

Ne  =  (2.2-4) 

whereupon  it  is  reflected  back  to  the  ground. 

At  Omega  frequencies  (-104  Hz),  the  critical  electron  density  (—1  electron/cm3)  given 
by  Eq.  2.2-4  occurs  at  an  altitude  of  50-60  km.  It  is  important  to  recognize  that  this  is  not  a 
sharply  defined  reflection  boundary.  A  closer  analysis  reveals  that,  at  the  altitude  of  the  critical 
electron  density,  the  wave  front  begins  to  turn  from  its  upward  direction,  propagate  horizontally, 
and  eventually  returns  to  earth.  Thus,  the  wave  actuals  scatters  through  the  ionosphere,  though 
it  may  appear  to  a  ground  observer  that  the  wave  is  r  iriy  reflected  from  a  sharply-defined 
surface  at  a  somewhat  higher  altitude  than  that  given  by  Eq.  2.2-4.  This  altitude  is  called  the 
effective  reflection  height  and  varies  from  about  70-75  km  (day,  sun  directly  overhead)  to  about 
85-90  km  (night,  at  antisolar  point). 

2.2.3  Photoionization/Chapman  Model 

In  descriptions  of  the  atmosphere/ionosphere,  a  length  scale  is  established  to  describe 
approximate  distances  over  which  various  processes  occur.  For  example,  if  the  atmosphere  is 
approximated  as  an  ideal  gas,  then  the  pressure  (P)  and  temperature  (T)  are  related  by 

P  =  nkT 


where  n  is  the  neutral  particle  number  density  and  k  is  Boltzmann’s  constant.  The  pressure  of  a 
column  of  gas  of  height  H  and  mass  density  g  in  the  (local)  earth’s  gravitational  field 
(described  by  the  gravitational  constant,  g)  is 

P  =  GgH 


where  g  =mn  and  m  is  the  average  mass  of  the  neutral  particles.  Equating  these  two  expressions 
for  the  pressure  gives  an  expression  for  H,  viz 


H 


kT 

mg 


5  H 


sc 


which  is  known  as  the  scale  height. 
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Solar  radiation  (photons)  striking  the  earth’s  atmosphere  loses  some  or  all  of  its  energy 
when  it  ionizes  neutral  particles.  This  absorption  of  solar  radiation  may  also  occur  without 
ionization  (particle  is  raised  to  an  excited  state).  A  measure  of  absorption/attenuation  of  solar 
radiation  is  the  optical  depth  (r )  which  is  given  by  (sun  directly  overhead) 

T  DoH^ 

where  n  is  neutral  particle  density  and  a  is  the  atomic  cross  section  for  absorption  (Ref.  7).  The 
optical  depth  is  generally  a  decreasing  function  of  altitude,  z  (following  n(z)),  and  the 
ionosphere  is  said  to  be  optically  thin  at  high  altitudes  and  optically  thick  at  low  altitudes.  It  is 
easily  shown  (Ref.  7)  that  if  IQ  is  the  radiation  flux  at  the  top  of  the  ionosphere,  then  the 
radiation  flux  at  altitude  z  is 

I(z)  =  Ioe'^ 

The  height  of  unit  optical  depth  is  the  height  (z0)  at  whi  :h  the  optical  depth  is  1,  i.e.  r  (z0)  =  1.  At 
this  altitude,  solar  radiation  flux  is  reduced  to  (l/'e)  of  its  incident  level.  This  height  represents  a 
balance  between  the  competing  processes  of  absorption  and  ionization:  At  lower  altitudes  more 
neutral  particles  are  available  for  the  solar  photons  to  ionize  (usually  into  ion-electron  pairs); 
however,  fewer  solar  photons  are  available  to  ionize  the  neutrals  due  to  absorption  of  solar 
photons  at  higher  altitudes. 

When  the  sun  is  not  directly  overhead  (x  >0°),  solar  radiation  traverses  more  of  the 
ionosphere  (thereby  intercepting  more  particles)  to  reach  the  same  altitude  and  the  optical 
depth  is  increased.  For  x  5  70°,  the  optical  depth  increases  by  a  factor  of  approximately  sec  %  ; 
for  larger  x  .  the  spherical  nature  of  the  atmospheric/ionosphere  must  be  accounted  for  and  the 
optical  depth  is  increased  by  the  Chapman  Function,  Ch  ( z,x  ,HSC)  (see  Appendix  A  for  details). 
Thus  the  optical  depth  is  more  generally  expressed  (for  a  fixed  ionizing  energy)  as 

r(z,  x)  =  n(z)cjHscCh(z,^,Hsc)  (2.2-5) 

The  electron  production  function,  introduced  in  Section  2.2.1  may  be  expressed  (for  a  fixed 
energy)  as 

q(z,  x)  =  —  oEoe-^)  (2.2-6) 

w 

where  w  is  the  energy  required  to  produce  one  ion-electron  pair  and  E0  is  the  energy  flux  at  the 
top  of  the  ionosphere,  (all  other  quantities  have  been  defined  earlier).  Because  e~r  increases 


and  n(z)  decreases  with  increasing  altitude,  the  electron  production  is  maximum  at  some 
intermediate  altitude. 

Based  on  the  definitions  and  assumptions  given  above,  Chapman  (Ref.  8)  showed  that, 
for  a  fixed  ionizing  energy  and  scale  height,  the  electron  production  function  can  be  written  as 

q(z’,X)  =  qmexpll-zVHsc-e-*'/^  Ch(z\z,Hsc)]  (2.2-7) 

where  qm  is  the  production  at  height  of  unit  optical  depth  and  z'  is  the  altitude  referenced  to  the 
height  of  unit  optical  depth. 

2.3  DAYTIME  PHASE  BEHAVIOR 

Based  on  characteristics  of  the  ionosphere  described  in  Section  2.2  and  the  principles  of 
wave  propagation,  a  relationship  is  established  in  this  section  between  signal  phase  variation 
and  solar  zenith  angle.  Included  in  the  relationship  are  necessary  (and  generally  well-known) 
ionospheric  and  propagation  parameters. 

2.3.1  Variation  of  Reflection  Height 

In  Section  2.2,  the  Chapman  electron  production  function  is  presented  (Eq.  2.2-7)  in 
terms  of  a  reduced  height  (relative  to  the  height  of  unit  optical  depth).  Eq.  2.2-7  may  be 
rewritten  as 

i 

loge(q/qm)  =  1  -  z'/Hsc  -  e'I'/HCh(z',x,  H*)  =  F(q(z',x)) 

The  function  F  is  a  monotonically  increasing  function  of  q;  thus  its  maximum  corresponds  to  q’s 
maximum.  The  maximum  F  (over  altitude)  is  easily  shown  to  occur  (ignoring  the  weak 
dependence  of  the  Chapman  function  on  z' )  when 

z'  =  Hsclogc  Ct^z'.z.Hsc)  =  z  — Zq  (2-3-1) 

where  z0  is  the  height  of  unit  optical  depth. 

For  those  situations  in  which  the  electron  density  changes  comparatively  slowly  over 
time,  dNe/dt  -  0  and  Eq.  2.2-3  becomes  (A  -  0) 


q  =  aeN^ 


In  this  case,  the  altitudes  for  maximum  electron  production  and  electron  density  coincide.  Even 
when  dNc/dt  cannot  be  neglected,  the  altitudes  of  the  corresponding  maxima  do  not  differ 
significantly  in  the  lower  ionosphere. 

Thus,  Eq.  2.3*1  essentially  describes  the  variation  of  the  peak  electron  density  altitude 
referenced  to  the  height  of  unit  optical  depth.  Although  the  effective  reflection  height  of  a  wave 
in  the  10-14  kHz  range  may  not  be  the  same  as  the  altitude  of  the  peak  electron  density,  the 
variation  of  the  two  defined  altitudes  with  solar  zenith  angle  is  expected  to  be  the  same  (Ref.  7). 
Hence  the  variation  in  reflection  height,  Ah ,  due  to  changing  solar  zenith  angle  behaves  as 

Ah  a  HsclogeChCz'.^Hsc)  (2.3-2) 

which  emphasizes  the  fact  that  the  variation,  Ah,  is  not  equal  to  the  variation  expressed  in 
Eq.  2.3-1  but  rather  that  behavior  with  changing  solar  zenith  angle  is  the  same. 

2.3.2  Phase  Change  Due  to  Variation  in  Reflection  Height 

Because  an  electromagnetic  wave  at  10  kHz  has  a  wavelength  nearly  equal  to  1/2  the 
distance  from  the  earth’s  surface  (where  the  wave  reflects)  to  the  ionospheric  region  where  the 
wave  scatters/retums  (-70  km  in  day),  the  signal  cannot  be  treated  as  a  “ray”  using  geometrical 
optics.  For  long-wavelength  signals,  two  descriptions  of  signal  propagation  are  generally  used, 
depending  upon  the  application.  The  wavehop  model  (Ref.  14)  treats  the  signal  received  at  a 
point  between  the  earth’s  surface  and  ionosphere  as  a  sum  of  wave  “hops”,  each  of  which  is  a 
distinct  integral  number  of  reflections  from  the  earth/ionosphere.  The  waveguide-mode  model 
(Ref.  15)  treats  the  earth's  surface  and  ionosphere  as  waveguide  boundaries  (the  “earth-iono¬ 
sphere  (El)  waveguide”)  and  the  signals  generated  and  propagated  inside  it  are  considered  to  be 
formed  from  a  series  of  “modes”  characteristic  of  the  waveguide.  The  latter  wave  propagation 
model  has  been  developed  extensively  (Ref.  12)  and  is  the  basis  for  most  of  the  semi-empirical 
modeling  described  in  this  report. 

Both  the  wavehop  and  the  waveguide-mode  models  lead  to  an  approximate,  linear 
relationship  between  signal  phase  variation  (A0)  on  signal  paths  within  the  El  waveguide  and 
variation  of  effective  reflection  height  (Ah)  with  solar  zenith  angle.  A  straightforward 
development  (Ref.  17)  yields 
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A 0  =  y  +  -^3)  Ah 

*  X  \  2Re  16Z3  I 


(2.3-3) 


where  L  is  the  path  length,  X  is  the  signal  wavelength,  z  is  nominal  reflection  height,  and  RE  is 
the  mean  earth  radius.  In  day,  z=70  km  so  that  (with  RE  =  6367  km) 

A0  =  y  (7.853  x  10-5  +  1.823  x  lO^Ah  (2.3-4) 

where  all  length  quantities  are  in  kilometers  and  A <p  in  cycles.  For  a  1  Mm  path  and  a  frequency 
of  10  kHz,  the  above  reduces  to 

A <p  =  8.084  x  10"3  Ah 

The  nominal  day-to-night  difference  in  effective  reflection  height  is  about  15  km  so  that 

A0day— night  —  12  cec/Mm  (10kHz) 

where  1  cec  =  0.01  cycle.  This  corresponds  approximately  to  the  diurnal  shift  observed  in  Ome¬ 
ga  phase  data  recordings  (Appendix  B). 

Substituting  Eq.  2.3-2  into  Eq.  2.3-3  yields 

A  0  =  CD  j(l  +  2.321  x  10'3A2)Hsc  l0gc  Ch(z\  x>  Hsc)  (2.3-5) 

where  all  lengths  are  in  kilometers  and  a  factor  of  7.853  x  10'5  is  taken  out  of  Eq.  2.3-4  and 
absorbed  into  the  unknown  constant  Cd  which  represents  the  proportionality  constant  implied 
by  Eq.  2.3-2.  Note  that  the  phase  variation  depends  upon 

•  frequency  (through  the  wavelength,  X ) 

•  solar  zenith  angle 

•  scale  height 

•  path  length. 

The  scale  height  used  for  most  atmospheric  models  (Ref.  9)  is  about  5-6  km  for  altitudes 
ranging  from  70  to  90  km.  A  default  of  Hsc  =  5  km  is  recommended  although  the  value  can  be 
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changed  if  so  demanded  by  the  observed  phase  data.  The  value  of  the  constant  Cd  is  to  be 
determined  by  the  data  as  discussed  in  Chapter  4.  The  functional  behavior  of  A <p  as  a  function 
of  x  given  by  Eq.  2.3-5  is  in  general  agreement  with  the  observed  data  (Appendix  B). 


2.4  NIGHTTIME  PHASE  BEHAVIOR 

Since  solar  radiation  is  not  directly  responsible  for  maintenance  of  the  ionosphere  at 
night,  the  development  in  this  section  does  not  strictly  follow  the  classical  Chapman  model 
described  in  Section  2.2.  The  basic  result  follows  from  semi-quantitative  arguments  and 
observational  data. 

2.4.1  Variation  of  Reflection  Height 

At  night,  direct  solar  photoionization  is  absent  and  the  ionospheric  profile  (electron 
density  vs  height)  changes  markedly.  The  D-region  effectively  disappears  at  night  primarily  due 
to  a  lack  of  strong  x-ray  and  extreme  ultraviolet  (EUV)  sources.  Thus  VLF  waves  are  scattered 
through  the  lower  E-region  (85-90  km)  which  is  maintained  by: 

•  Cosmic  rays 

•  Galactic  radiation  (x-ray  stars,  etc) 

•  Scattered  solar  radiation. 

Of  these  three  sources,  the  most  important  is  the  scattered  solar  radiation  which  consists  of 
three  principal  types  (Ref.  9):  (1)  Lyman-a,  (2)  Lyman-/?,  and  (3)  He  H.  In  the  lower  E-region, 
Lyman-a  fluxes  predominate  over  the  other  sources.  Lyman-a  (the  primary  transition  between 
the  first  two  orbitals  of  hydrogen)  is  produced  in  immense  quantities  by  the  sun  but  it  has  a 
relatively  long  wavelength  (1216  A)/low  energy  so  it  does  not  normally  penetrate  the  lower 
ionosphere  even  on  the  dayside.  However,  nitric  oxide  in  the  ionosphere  is  very  efficiently 
ionized  by  Lyman-a  photons  so  that  even  small  fluxes  contribute  to  D/E-region  ionization. 
Lyman-a  radiation  also  passes  through  the  earth’s  geocorona  and  resonantly  scatters  from  the 
large  local  population  of  hydrogen  atoms/ions.  This  scattering  is  relatively  efficient  and  can 
produce  significant  fluxes  of  Lyman-a  radiation  incident  on  the  nightside  ionosphere  (Ref.  9). 
Satellite  measurements  of  this  radiation  (Ref.  10)  show  a  distinct  solar  zenith  angle  dependence 
(0  <  x  <  n)  Figure  2.4-1  shows  the  Lyman-a  radiation  received  at  a  near- zenith  angle  from 
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SOLAR  ZENITH  ANGLE  (deg) 

Figure  2.4-1  Theoretical  and  Experimental  Near-zenith  Lyman-alpha  Intensity 
as  a  Function  of  Solar  Zenith  Angle  (from  Ref.  10) 

the  satellite.  This  behavior,  together  with  the  nadir  intensity  behavior  (see  Ref.  10)  indicates  an 
approximate  cosine-type  functional  dependence  on  solar  zenith  angle. 

If  the  Lyman-a  radiation  actually  controls  the  nighttime  effective  reflection  height,  then 
the  observed  signal  phase  should  reflect  this  behavior  during  the  nighttime  measurement  period 
(phase  should  peak  when  the  average  solar  zenith  angle  on  a  path  is  maximum).  If  no  effective 
source  of  nighttime  ionization  is  present  but  the  recombination  processes  are  very  slow,  the 
phase  should  gradually  rise  during  the  time  the  ionosphere  controlling  the  signal  path  is  fully 
dark.  The  observed  data  (Appendix  B)  frequently  shows  a  slight  peak  during  the  nighttime 
portions  indicating  a  solar  zenith  angle  dependence. 

In  order  to  derive  a  quantitative  relationship  between  Lyman-a  flux  and  solar  zenith 
angle,  first  assume  that  the  effective  flux  is  incident  normal  to  the  nignttime  ionosphere  (i.e., 
radially  inward  at  any  point  on  the  nightside).  This  corresponds  to  %  =  0  in  the  Chapman 
formalism  and  Eq.  2.2- /  may  be  written 

q(z',0)  =  qm  exp[l  -  z'/Hsc- e‘z /»*]  (2.4-1) 
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since  Ch(z,0,Hsc)  =  sec(O)  =  1.  As  noted  earlier,  qm  is  the  electron  production  relative  to  the 
height  of  unit  optical  depth  (zo),i.e.,  qm  =  q(z'  =  0,0)  where  z'  =  z-zo  .  Using  Eqs.  2.2-5 
and  2.2-6,  qm  may  be  written  (x  =  0) 

TE0e~r|  _  E0 
wHv  I  i  ewH 


Qm  — 


sc 


lsc  lT  =  1 

Since  w  is  the  minimum  energy  required  to  create  an  ion-electron  pair,  E0/w  is  just  the  ionizing 
efficiency  ( rj )  times  the  incident  photon  flux  (I„)  at  the  top  of  the  ionosphere.  Thus 


Qm  — 


eH 


(2.4-2) 


SC 


This  differs  from  the  solar  Chapman  model  in  that  the  radiation  is  normally  incident  but  the 
external  flux  (or  intensity!  IQ  is  a  function  of  solar  zenith  angle. 

Further  assume  that  the  ionosphere  is  approximately  in  equilibrium  so  that  Eq.  2.2-3 
becomes  (with  A  =  0) 


q  =  «eN? 


which,  when  combined  with  Eqs.  2.4-1,  2.4-2  yields 

«,N§  =  ^  expf-z'/H*-e-’VH.l  (2.4-3) 

Hsc  L  J 

The  reasonable  assumption  is  now  made  that  wave  reflection  occurs  at  a  fixed  value  of 
electron  density,  Ne  (i.e.,  a  fixed  plasma  frequency;  see  Eq.  2.2-4).  Thus,  in  Eq.  2.4-3,  only  the 
reduced  height  z'  changes  in  response  to  the  solar  zenith  angle-dependent  intensity  ID.  Since  the 
height  of  unit  optical  depth  for  Lyman-a  radiation  is  75  km  (Ref.  9)  and  the  nighttime  reflection 
height  is  80-90  km, 

1  <  z'/Hsc  <  3 

and  the  exponential  term  inside  the  brackets  in  Eq.  2.4-3  may  reasonably  be  neglected  in 
comparison  to  the  linear  term.  Thus,  Eq.  2.4-3  becomes 

e'V-a.N?  = 


or 

z'  =  Hsc  logeloCt)  +  (terms  independent  of  *)  (2.4-4) 

Recall  that  the  satellite  observations  (see  Fig.  2.4-1)  showed  that  the  logarithm  of  Lyman-a 
radiation  intensity  has  a  cosine-type  dependence  on  solar  zenith  angle.  This  fact,  coupled  with 
Eq.  2.4-4  implies  that  z'  =  z-Zq  depends  linearly  on  cos*  and  hence  it  follows  that 

Ah  =  A  cos*  +  B  (2.4-5) 

where  it  is  assumed,  that  the  variation  in  altitude  of  the  ionospheric  level  having  constant 
electron  density  (critical  value  for  reflection  of  10-14  kHz  waves)  referenced  to  the  height  of 
unit  optical  depth  is  the  same  as  the  variation  in  reflection  height.  In  Eq.  2.4-5,  A  and  B  are 
constants  to  be  determined  from  the  observed  phase  data  and  by  matching  phase  at  the  day/ 
night  boundary. 

2.4.2  Phase  Change  Due  to  Variation  in  Reflection  Height 

The  signal  phase  variation  due  to  the  reflection  height  change  given  by  Eq.  2.4-5  is 
obtained  from  the  general  form,  Eq.  2.3-3.  Thus,  with  z  =  87  km  for  the  nominal  nighttime 
ionospheric  reflection  height, 

A#  =  j  (1  +  1.209  x  10“3A2)(A  cos*  +  B)  (2.4-6) 

Note  that  the  coefficient  A  must  be  negative  since  the  phase  increases  as  the  reflection  height 
increases  and  thus  as  the  solar  zenith  angle  increases  from  n/2  to  n  . 


2.5  ONSET/RECOVERY  TIME  PARAMETERS 

The  two  previous  sections  addressed  phase  behavior  as  a  function  of  solar  zenith  angle 
for  local  daytime  and  nighttime  ionospheric  conditions.  Accurate  specification  of  the  solar 
zenith  angle  “boundary”  between  these  two  temporal  domains  is  important  for  phase  prediction 
accuracy  as  illustrated  in  Fig.  2.1-2.  Clearly,  a  scheme  which  predicts  the  correct  shape  of  a 
diurnal  phase  trace  but  is  displaced  by  20-30  cec  due  to  an  onset  time  prediction  error  is  not  a 
useful  predictor. 
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2.5.1  Definition  of  Onset/Recovery  Time  Parameters 

Sunrise/sunset  onset/recovery  times  refer  to  the  phase  behavior  along  an  entire  signal 
path  defined  by  a  transmitting  source  at  one  end  and  a  receiver  at  the  other  end.  Since  results 
from  the  waveguide-mode  model  are  used  for  the  spatial  sub-models,  the  signal  path  along  the 
two-dimensional  El  waveguide  is  broken  into  path  segments  (see  Fig.  2.5-1).  Within  each  of 
these  0.01-radian  segments,  the  path  spatial  properties  are  assumed  constant  and  the  waveguide 
is  homogeneous.  Thus,  local  ionospheric  properties  are  defined  at  the  segment  level  and  path 
characteristics  reflect  the  collective  segment  level  properties. 

A  local,  daytime  ionosphere  is  defined  as  one  whose  ionization  processes  are  solar- 
controlled.  Clearly,  solar  radiation  controls  the  ionization  of  the  D-region  for  0  <  %  <  Jt/2.  For  % 
>  nil,  solar  radiation  must  pass  through  lower  altitudes  and  re-emerge  at  the  desired  altitude 
(see  Fig.  2.5-2).  This  property,  which  is  embodied  in  the  appropriate  approximations  to  the 
Chapman  function,  means  that  the  rate  of  attenuation  of  solar  radiation  (due  to  changing  solar 
zenith  angle)  becomes  very  rapid.  Solar  photoionization  is  still  expected  to  dominate  other 
sources  until  the  radiation  “beam”  incident  upon  the  75-85  km  ionospheric  layer  (from  below) 


Figure  2.5-1  Segmentation  of  the  Signal  Path 
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attains  an  enroute  minimum  of  about  50  km.  At  this  point  the  concentration  of  ozone  which 
efficiently  absorbs  EUV  radiation  becomes  substantial  and  ionization  due  to  solar  radiation 
effectively  ceases.  The  solar  zenith  angle  corresponding  to  this  condition  is  called  the  threshold 
solar  zenith  angle  CtTHR )  and,  for  a  nighttime  reflection  height  of  about  85  km,  is 
approximately  96°  (see  Fig.  2.5-2). 

Ionization  does  not  occur  immediately  after  radiation  is  incident  upon  a  region  of  neutral 
constituents.  Solutions  of  the  continuity  equation,  Eq.  2.2-3,  for  an  impulse-type  electron 
production  function  show  that  the  electron  density  “impulse”  follows  the  impulsive  production 
by  a  time  interval,  r,  given  by 

1  1 

r  = - =  — t=  —  lonosphenc  response  time 

2aE'Nc  2j^ 

In  the  E-region  the  recombination  coefficient  is  given  as  10~7  cm3/sec  (Ref.  13).  When  used  with 
D-region  electron  densities,  however,  this  numerical  value  of  aE  yields  an  unacceptably  large 


SOLAR 

RADIATION 


Ozone  layer  cutoff  =  50  km 

Reflection  height  2  85  km 

Threshold  solar  zenith  angle  (Xthr)  =  96° 

Figure  2.5-2  Illustration  of  the  Threshold  Solar  Zenith  Angle  which  Separates 
Day  and  Night  for  the  Lower  Ionosphere 
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value  of  r  .  Observational  data  (Ref.  11)  suggests  ionospheric  response  times  of  about  1/8  hour 
(day)  and  1.13  hours  (night).  Fig.  2.5-3  illustrates  the  determination  of  ionospheric  response 
time. 

In  determining  the  onset/recovery  of  daytime  or  nighttime  conditions  along  a  signal  path, 
the  interaction  of  the  wave  with  the  ionosphere  must  be  considered  in  addition  to  the 
ionospheric  changes.  Since  the  wave  interacts  with  the  ionosphere  over  a  substantial  fraction  cf 
the  path,  it  is  important  to  know  the  point  on  the  path  at  which  the  first  effective  path/ionosphere 
interaction  takes  place. 

An  Omega  antenna  can  be  modeled  as  monopole  source  which  radiates  energy  outward 
and  upward  from  its  location  on  the  earth’s  surface.  The  steeply  launched  portions  of  the 
transmitted  wave  interact  first  with  the  ionosphere  but  are  reflected  between  the  earth’s  surface 
and  ionosphere  before  traveling  outward  to  any  appreciable  distance.  Those  portions  of  the 


D-Region  Balance  Equation 

dN« 


df 


=  q-aNe 


Ne  =  Electron  density 
q  =  Electron  production 
a  =  Recombination  coefficient 


Early  Experimental 
Results  Indicated 

x0  =  6  minutes  (day) 

tn  =  68  minutes  (night) 


Figure  2.5-3  Origin  of  the  Day/Night  Ionospheric  Response  Time 
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transmitted  wave  which  are  launched  nearly  tangentially  to  the  earth's  surface  travel  farthest 
outward  before  interacting  with  the  ionosphere.  Because  of  its  relatively  low  attenuation,  this 
wave  structure  generally  dominates  all  others  (in  signal  energy)  at  distances  (along  the  earth’s 
surface)  beyond  the  wave’s  first  interaction  with  the  ionosphere.  The  wave  structure  which 
survives  in  the  far  field  is  launched  at  an  angle  slightly  less  than  90°  from  the  zenith  which 
satisfies  the  reinforcement  condition  between  the  reflected  wave  and  the  direct  wave. 
Satisfaction  of  this  condition  insures  that  the  wave  is  an  eigenmode  (Mode  1)  of  signal 
propagation  (at  a  particular  Omega  frequency)  in  the  earth-ionosphere  waveguide.  Further 
discussion  of  these  modes  is  given  in  Chapter  3. 

Figure  2.5-4  illustrates  the  geometry  of  the  Mode  1  wave  interaction  with  the  ionosphere. 
The  point  of  first  ionospheric  interaction,  P,  is  separated  from  the  transmitting  source  by  an 
internal  angle,  a ,  measured  at  the  center  of  the  earth.  Since  Re«  >  >  h,  VLF  signal  paths  are 
usually  measured  by  their  arc  length  over  the  earth’s  surface.  The  table  in  Fig.  2.5-4  indicates 
how  a  varies  with  ionosphere  illumination  and  frequency.  Night  values  are  larger  than  day 
because  h  is  greater  so  that  the  wave  travels  outward  further  before  intersecting  the  ionosphere. 
Angular  values  are  higher  for  13.6  kHz  than  10.2  kHz  because  the  corresponding  eigenangle  is 
larger  for  the  higher  frequency  (wavelength  is  smaller  so  mode  resonance  condition  more 

o 

closely  approximates  a  rp  -  90  launch).  Because  of  the  general  reciprocity  features  of 
electromagnetic  waves  (reciprocity  is  not  strictly  valid  due  to  the  anisotropic  magnetized 
ionospheric  plasma),  the  same  geometry  holds  for  the  receiving  end  of  the  path.  Thus,  if  the 
transmitter  in  Fig.  2.5-4  is  replaced  by  the  receiver,  then  P  becomes  the  last  point  of  ionospheric 
interaction  for  the  Mode  1  signal  (angular  relationships  are  approximately  the  same). 

2.5.2  Numerical  Testing  of  Sub-model  Forms 

Since  the  diurnal  sub-model,  by  definition,  contains  no  provision  for  any  sources  of 
phase  variation  other  than  the  diurnal  variation  of  effective  ionospheric  reflection  height, 
comparison  with  observed  data  on  paths  of  conventional  length  is  not  feasible,  primarily 
because  of  the  coupling  between  the  diurnal  and  spatial  sub-models.  However,  the  onset/ 
recovery  time  parameters  are  essentially  independent  of  the  spatial  contributions  to  the  phase 
prediction  and  thus  can  be  tested. 

Day  and  night  functional  forms  presented  in  Sections  2.3  and  2.4  are  computed  at  each 
segment  (depending  on  solar  zenith  angle)  and  combined  with  the  reference  phase  depending 
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p 

V  -  Real  part  of  Mode  1  eigenangle  foi 
frequency /illumination  condition 

h  -  Ionospheric  reflection  height 

Re  =  Radius  of  earth 

a  *  Angular  extent  of  near-field  region 

(no  ionospheric  interaction  with  Mode  1  signal) 

P  =  Point  cf  first  ionospheric 
interaction 

a  (radians) : 


Figure  2.5-4  Schematic  of  Path  Geometry  to  Illustrate  Point  of  First  Ionospheric 
Interaction  for  a  Mode  1  Signal 

only  on  path  length.  The  path/time  computer  program  (developed  under  this  effort)  contains 
path  segmentation  and  solar  position  subroutines  similar  to  those  used  in  earlier  PPC  models 
(Ref.  3).  Required  inputs  to  the  program  are  transmitting  station,  receiving  site,  signal 
frequency,  and  time  (month/day/year).  Input  values  are  also  used  to  identify  the  path  segment 
number  defining  point  of  first/last  ionospheric  interaction,  solar  zenith  angle  threshold,  and 
day/night  ionospheric  response  times.  The  program  computes  path  phase  values  at  six-minute 
intervals  based  only  on  diurnal  effects  and  reference  wave  number.  Figure  2.5-5  shows  a 
sample  printout  from  the  program  which  lists  the  beginnings  of  sunset,  night,  sunrise,  and  day 
(note  that  path  sunset  recovery  time  is  the  same  as  the  beginning  of  path  night  and  path  sunrise 
recovery  time  is  the  same  as  the  beginning  of  path  day).  Time  units  are  in  0.1  hour  (6  minutes) 
and  the  “-1”,  “-2”  for  each  time  domain  allows  for  double-diurnal  effects  (see  Appendix  B). 
The  phase  data  given  at  the  bottom  of  the  printout  for  each  6  minute  interval  is  plotted  in 
Fig.  2.5-6. 
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****  THIS  IS  THE  OATA  FOR  PLOT  FILE  BA260988.13T  **** 
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Figure  2.5-5  Sample  Printout  from  Path/Time  Computer  Program  to  Test 
Diurnal  Sub-model  Forms 
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Figure  2.5-6  Plot  of  Diurnal  Phase  Prediction 

It  is  important  to  recognize  that  the  diurnal  phase  behavior  plotted  in  Fig.  2.5-6  does  not 
account  for  any  spatial  variation  in  phase  (i.e.,  that  due  to  the  geomagnetic  field  or  ground 
conductivity)  and  thus  cannot  be  properly  compared  to  observed  data.  It  is  shown  to  indicate 
how  phase  variation  would  appear  in  a  spatially  homogeneous  world. 

Comparison  of  predicted  and  observed  onset/recovery  times  for  sunrise  and  sunset  is 
given  in  Appendix  B.  There  it  is  shown  that  the  predictions  are  in  close  agreement  with  the 
observations  studied  using  default  onset/recovery  time  parameters. 

2.6  SUMMARY 

In  this  chapter  the  specific  functional  behavior  of  the  VLF  signal  phase  with  solar  zenith 
angle  is  developed  for  two  broad  classes  of  local  ionospheric  conditions: 

•  control  by  direct  solar  radiation  (day) 

•  control  by  nighttime  radiation  sources  (night). 

The  daytime  phase  variation  is  based  on  a  Chapman  model  of  electron  photoionization/absorp¬ 
tion  and  numerical  values  indicate  a  qualitative  agreement  with  observed  data.  The  nighttime 
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phase  variation  is  developed  from  empirical  measurements  of  scattered  Lyman-alpha  flux  using 
a  method  of  analysis  similar  to  Chapman’s.  The  resulting  phase  variation  exhibits  a  phase 
change  with  solar  zenith  angle  at  night  —  a  phenomenon  also  detected  in  the  observed  data. 

Parameters  defining  the  path  sunrise/sunset  onset/recovery  times  are  identified  and  a 
path/time  computer  program  is  written  to  compare  the  predicted  times  with  observed  data. 
Comparison  of  the  predicted  times  using  default  parameters  with  a  limited  sample  of 
observations  show  excellent  agreement.  Although  eventual  calibration  with  observed  phase 
data  will  give  more  precise  values,  these  results  show  that  the  default  values  are  good  estimates 
of  the  onset/recovery  time  parameters. 
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3. 


GEOMAGNETIC/CONDUCTIVITY  SUB  MODEL 


3.1  INTRODUCTION 

This  chapter  presents  the  methodology  followed  to  develop  the  structure  of  the  geomag¬ 
netic/conductivity  sub-model.  This  sub-model  is  the  spatial  component  of  the  overall  Omega 
signal  phase  prediction  model  (referred  to  as  the  1990  PPC  model)  described  in  this  report.  The 
1990  PPC  model  predicts  the  phase  of  a  received  Omega  signal  (a  multimode  signal)  by  assum¬ 
ing  the  signal’s  Mode  1  component  amplitude  greatly  exceeds  all  other  component  mode  ampli¬ 
tudes. 

The  Mode  1  signal  phase  can  be  separated  into  two  components: 

0  =  0l(d)  +  02 

where  0i  depends  on  path  length  (d)  and  fa  does  not.  The  first  component  may  be  expressed  as 
(assuming  insignificant  signal  dispersion): 


where  k  is  the  wave  number,  f  is  the  frequency  and  v  is  the  phase  velocity.  For  a  signal  of  fixed 
frequency  and  signal  path  of  fixed  length,  fa  thus  depends  on  wave  number  or  phase  velocity. 
For  this  reason,  fa  is  called  the  wave  number  or  phase  velocity  component.  The  component  fa  is 
called  the  excitation  factor  phase.  This  sub-model  describes  the  theoretical  dependence  of  these 
Mode  1  signal  phase  components  as  functions  of  the  locally-varying  spatial  parameters  (i.e., 
ground  conductivity,  geomagnetic  field  in  the  ionosphere,  and  signal  propagation  direction)  at 
each  segment  of  the  signal  path  for. 

•  Each  of  the  two  extreme  (day  and  night)  solar  illumination  conditions 

•  Each  of  three  (10.2,  11  1/3,  and  13.6  kHz)  Omega  navigational  frequen¬ 
cies. 

An  overview  of  the  Omega  signal  propagation  mechanism/theory  is  given  in  Section  3.2. 
The  methodology  used  for  developing  the  theoretical  database,  from  which  the  sub-model  is 
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derived,  is  given  in  Section  3.3.  The  methodology  used  to  develop  the  sub-model  and  the  result¬ 
ing  sub-model  structure  is  presented  in  Section  3.4.  The  chapter  is  summarized  ir  Section  3.5. 


3.2  OMEGA  SIGNAL  PROPAGATION  MECHANISM/THEORY 

As  presented  in  Chapter  2,  Omega  signals  propagate  in  the  El  waveguide,  i.e.,  the  space 
between  the  earth’s  surface  and  the  D-region  of  the  ionosphere.  In  this  waveguide,  signals  ra¬ 
diated  by  a  transmitting  station  propagate  essentially  along  great-circle  radials  emanating  from 
the  station.  The  propagation  characteristics  of  an  Omega  signal  along  a  radial  path*  are  func¬ 
tions  of  the  electromagnetic  properties  of  the  path  waveguide.  The  lower  boundary  of  this  wave¬ 
guide  is  the  earth’s  surface  characterized  by  a  spatially-varying  ground  conductivity  (and 
associated  permittivity).  Table  3.2-1  lists  the  worldwide  values  of  the  earth’s  ground  conductiv¬ 
ity  (and  permittivity)  at  VLF/Omega  signal  frequencies. 

The  upper  boundary  of  the  El  waveguide  is  the  D-region  of  the  ionosphere  which  is  ani¬ 
sotropic  due  to  the  vector  geomagnetic  field.  Most  VLF  propagation  prediction  models  (in  par¬ 
ticular,  the  waveguide-mode  model  (Ref.20))  assume  the  electron  density  varies  exponentially 
with  altitude  in  the  vicinity  of  the  reflection  height,  i.e., 

Ne(z)  =  Ne(zo)  e  b*z“z°) 

Table  3.2-1  Earth’s  Ground  Conductivity  Levels  and  Associated 
Conductivity/Permittivity  Values  (Ref.  19) 


LEVEL 

CONDUCTIVITY 

(mho/m) 

RELATIVE 

PERMITTIVITY 

1  (Ice  Cap) 

1.0  X  10'5 

10.0 

2 

3.2  x  10'5 

10.0 

3 

1.0  x  10'4 

15.0 

4 

3.2  x  10'4 

15.0 

5 

1.0  x  10'3 

15.0 

6 

3.2  x  10'3 

15.0 

7 

1.0  x  10'2 

20.0 

8 

3.2  x  10-2 

20.0 

9 

1.0  x  10’1 

45.0 

10  (Seawater) 

4.0 

81.0 

•“Path"  refers  to  a  two-dimensional  waveguide  in  which  one  dimension  is  the  altitude 
above  the  earth’s  surface  and  the  other  dimension  is  the  arc  length  from  the  transmitting 
station  along  the  surface  of  the  earth. 


where  z<j  is  the  effective  reflection  height  and  b  (b>0)  is  the  electron  density  gradient.  These  mod¬ 
els  also  assume  that  the  collision  frequency  profile  is  exponential  in  the  vicinity  of  the  effective 
reflection  height,  i.e., 

ve(z)  =  ve(zo)  e  _a(z"Zo> 


where  ve(z)  is  the  electron  collision  frequency  (electron/neutral  and  electron/ion)  and  a  (a>0)  is 
the  collision  frequency  gradient.  For  very  low  signal  frequencies,  the  direct  conductivity  is  de¬ 
fined  as 


oQ(z)  = 


QJp(z) 

ve(z) 


where  cop  is  the  plasma  frequency  defined  in  Chapter  2. 

In  Section  2.2  it  is  shown  that  ouip  is  proportional  to  the  electron  density,  Ne(z).  Thus  the 
direct  conductivity  is  proportional  to  Ne(z)/ve(z)  and  hence  is  also  exponential  (as  a  function  of 
altitude)  in  the  vicinity  of  z  -  Zq  with  gradient  fi  where 

=  a  +  b 

The  effective  reflection  height,  Zq  ,  and  the  conductivity  gradient,  p ,  are  the  parameters  re¬ 
quired  by  waveguide-mode  model  calculations  of  signal  propagation.  Numerical  values  of  these 
parameters  are  given  in  Table  3.2-2. 


Since  Omega  signal  path  parameters  vary  both  spatially  (a.ung  the  path)  and  temporally 
(at  each  point  along  a  path),  Omega  signal  paths,  and  hence  the  associated  waveguides,  are  in¬ 
homogeneous.  The  waveguide-mode  model  transforms  an  inhomogeneous  path  waveguide  into 
a  set  of  homogeneous  waveguides  by  means  of  the  following  procedure:  (1)  divide  the  inhomo¬ 
geneous  waveguide/path  into  a  large  number  of  concatenated  homogeneous-property  segments, 
(2)  treat  each  segment  as  a  homogeneous,  infinitely  long  waveguide  with  properties  correspond¬ 
ing  to  the  upper  and  lower  boundaries  of  the  segment,  and  determine  the  allowed  modes/signal 


Table  3.2-2  Typical  Values  of  Daytime  and  Nighttime  Ionospheric  Parameters 


IONOSPHERE 

ILLUMINATION 

CONDITION 


Day 

Night 


IONOSPHERIC 
REFLECTION  HEIGHT 
(km) 


70 


IONOSPHERIC 
CONDUCTIVITY  GRADIENT 
(km*1) 


0.3 


87 


0.5 


behavior,  (3)  determine  mode  conversion  effects  (scattering  of  the  incident  signal)  occurring  at 
each  of  the  interfaces  between  the  homogeneous  segments,  and  (4)  connect  the  multimode  sig¬ 
nal  solutions  at  each  segment  by  properly  accounting  for  the  mode  conversion  effects  at  the  in¬ 
terfaces. 

In  the  waveguide-mode  theory,  an  electromagnetic  signal  field  is  expressed  as  a  sum  of 
the  signal  fields  corresponding  to  the  “characteristic”  modes  of  the  waveguide.  The  modes  char¬ 
acterize  the  solutions  of  the  waveguide-mode  equation,  which  describes  the  allowable  wave 
launch  angles  that  establish  reinforcement  between  the  transmitted  and  reflected  waves.  Each 
mode  has  an  associated  complex  eigenvalue  (also  referred  to  as  an  eigenangle).  The  eigenangle 
is  a  function  of  the  path  parameters  and  signal  frequency.  The  signal  parameters  of  a  mode,  i.e., 
the  amplitude  and  phase,  are  functions  of  the  following  components: 

•  Attenuation  rate  —  the  rate  at  which  the  amplitude  of  the  signal  changes 
along  the  path 

•  Phase  velocity  —  the  rate  at  which  a  point  of  constant  phase  appears  to 
propagate  along  the  path 

•  Excitation  factor  (a  complex  quantity)  —  the  relative  efficiency  with 
which  the  signal  is  excited  by  the  transmitting  antenna,  or  received  by 
the  receiving  antenna. 

These  signal  components  are  functions  of  the  mode  eigenvalue  and  signal  frequency.  The  atten¬ 
uation  rate  and  excitation-factor  amplitude  jointly  determine  the  amplitude  of  the  mode  signal. 
The  phase  velocity  and  excitation-factor  phase  jointly  determine  the  phase  of  the  mode  signal. 

Individual  mode  solutions  of  the  waveguide-mode  equation  are  distinguished  from  each 
other  by  numbering  with  positive  integers.  The  lowest  phase-velocity  mode  is  usually  labeled  as 
“Mode  1”;  Mode  2  has  the  next  lowest  phase  velocity,  and  so  on.  In  this  mode  numbering  sys¬ 
tem,  odd-  and  even-numbered  modes  are  the  transverse-magnetic  (TM)  and  transverse-electric 
(TE)  modes,  respectively.  The  electric  field  of  a  TM  (TE)  mode  is  vertically  (horizontally)  polar¬ 
ized.*  Therefore,  the  TM  mode  signals  are  most  effectively  received  by  antennas  such  as  a  verti¬ 
cal  whip  or  a  loop  antenna,  whose  symmetry  axis  is  parallel  to  the  ground. 

Mode  1  signals,  in  addition  to  having  the  lowest  phase  velocity,  usually  have  the  smallest 
attenuation  rate  and  a  relatively  high  excitation-factor  amplitude.  As  a  consequence,  Omega 
station  signals  are  usually  dominated  by  their  respective  Mode  1  components.  Omega  signals 

'Vertical,  i.e.,  perpendicular  to  the  earth’s  surface. 
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are  not  dominated  by  their  associated  Mode  1  components  in  the  station  “near-field”  region 
which  typically  extends  to  ranges  of  300-500  km  during  daytime,  1000-3000  km  during  night¬ 
time,  and  for  all  ranges  along  certain  westerly-directed,  nighttime  paths.  The  signal  condition  in 
which  Mode  1  is  not  dominant  is  commonly  referred  to  as  the  “modal  condition.”  A  modal  sig¬ 
nal  is  either  composed  of  several  semi-dominant  modes  or  dominated  by  a  single  mode,  other 
than  Mode  1.  A  modal  signal  is  unsuited  for  Omega  navigation. 


3.3  THEORETICAL  DATABASE 


This  section  presents  the  approach  used  to  generate  the  theoretical  database  required  to 
develop  the  geomagnetic/conductivity  sub-model.  Also  described  in  this  section  are  relevant 
features  of  the  Mode  1  database  which  consists  of  theoretical  values  of  the  Mode  1  signal  phase 
components  (i.e.,  phase  velocity  and  excitation-factor-phase)  as  functions  of  the  worldwide  val¬ 
ues  of  the  spatial  parameters  (i.e.,  ground  conductivity  and  geomagnetic  field).  The  Mode  1  sig¬ 
nal  phase  components  are  computed  for: 

•  Daytime  and  nighttime 

•  The  three  signal  frequencies:  10.2,  11  1/3,  and  13.6  kHz. 


The  conductivity  at  a  path  segment  is  simply  the  ground  conductivity  (read  from  a  conductivity 
map/file)  at  the  center  of  the  path  segment.  Other  spatial  parameters  include  the  vector  geo¬ 
magnetic  field  and  the  direction  of  signal  propagation.  If  the  geomagnetic  field  is  approximated 
as  an  earth-centered  magnetic  dipole  field,  then,  in  terms  of  path  coordinates  (i.e.,  two  axes  par¬ 
allel  (T||)  and  perpendicular  (Ti)  to  the  path  in  the  local  tangent  plane  of  the  earth,  and  one 
axis  radially  upward),  the  field  components  may  be  written 

Bja  =  Msin^mcos^m  ;  BtI(  =  Mcos/?mcos0m  ;  Br  =  2Msin0m 


where 


M - 'Be°' 

(1  +  3  sin2  0m)1/7 


—  geomagnetic  path  bearing  ;  9m  =  geomagnetic  latitude 
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Here,  |  Beq  |  is  the  magnitude  of  the  geomagnetic  field  at  the  earth’s  surface*  at  the  geomagnet¬ 
ic  equator.  Thus  only  and  (in  addition  to  the  fixed  constant  j  Beq  | )  are  needed  to  specify 
the  geomagnetic  field  components  relative  to  the  path  orientation. 

As  described  in  Section  3.2,  the  phase  components  of  a  mode  signal  are  derived  from  the 
eigenvalue  of  the  mode  and  signal  frequency.  Thus  to  form  the  required  database  for  the  sub¬ 
model,  a  database  of  Mode  1  eigenvalues  must  be  computed  which  is  parametric  in  the  spatial 
path  parameters,  signal  frequency,  and  illumination  condition.  The  approach  used  here  to  gen¬ 
erate  such  an  eigenvalue  database  is  to  use  the  IP P  computer  code  (Ref.  18),  together  with  the 
MODESRCH  computer  code  (Ref.  19)  both  developed  at  the  Naval  Ocean  Systems  Center. 

Since  the  desired  database  currently  exists  for  10.2  and  13.6  kHz  (Ref.  21),  the  11-1/3 
kHz  database  was  generated  under  this  effort.  The  11-1/3  kHz  database  was  generated  using 
IPP,  aided  by  MODESRCH  (as  needed)  for  conductivity  level  10  only;  otherwise  the  database  in¬ 
cluded  values  for  the  same  path  parameter  combinations  as  the  10.2  and  13.6  kHz  databases. 
The  combinations  included  each  of  the  ten  ground  conductivity  levels,  at,  (for  10.2  and 
13.6kHz  only),  the  geomagnetic  latitude,  9m,  and  geomagnetic  bearing,  /5m,  as  listed  in 
Table  3.3-1.  The  geomagnetic  latitude/bearing  grid  used  in  the  nighttime  database  computa¬ 
tions  is  purposely  chosen  to  be  non-uniform  with  increasing  number  of  grid  computation  points 


Table  3.3-1  Geomagnetic  Latitude  and  Bearing  Grids  used  in  the  Database 


SPACING 

GRID  BOUNDS 

ILLUMINATION 

GRID 

TYPE 

BETWEEN  GRID 
LINES  A  dm  X  A£m 

(deg) 

GEOMAGNETIC 
LATITUDE,  0m 
(deg) 

GEOMAGNETIC 

BEARING. 

(deg) 

Degree 

D-l 

10  x  10 

0  £  9m  <  90 

90  £  /3m  <;  270 

Night 

N-l 

10  x  10 

30  £  9m  £  90 

90  s  /Sm  <;  270 

N-2 

5x5 

10  <;  9m  <  30 

90  <;  /9m  ^  270 

N-3 

1  x  5 

0  ^  9m  <  10 

90  <;  pm  <,  170 

and 

210  <;  /?m  s  270 

N-4 

1  x  1 

0  £  9m  <  10 

170  <  <  210 

•There  is  negligible  difference  between  |B£Q|at  the  earth’s  surface  and  |Beq| 
at  D-region  altitudes. 


at  lower  geomagnetic  latitudes  and  westerly  geomagnetic  bearings  where  mode  eigenvalues 
change  very  rapidly  for  a  small  change  in  the  spatial  parameters. 

The  available  10.2  and  13.6  kHz  databases  include  three  daytime  modes  (Modes  1,  2,  3) 
and  five  nighttime  modes  (Modes  1,  2,  3,  4,  X).  The  11-1/3  kHz  database  is  limited  to  two 
daytime  modes  (Modes  1, 2)  and  three  nighttime  modes  (Modes  1, 2,  X).  Because  of  the  approx¬ 
imate  symmetry  of  the  eigenvalues  (e)  i.e., 

e  (0m, An)  =  e  Mm, An)  (3-3-1) 

tfanM  =  e  (0m,  180° -/?m)  (3.3-2) 

the  eigenvalues  were  computed  (and  are  contained  in  the  database)  only  for  0°  <  9m  S  90°, 
and  90°  <  <  270°. 

Figures  3.3-l(a)  and  3.3-l(b)  depict  the  9m-  and  /?m -dependence  of  the  eigenvalues  for 
daytime 'nighttime,  13.6  kHz  signal  frequency,  and  conductivity  level  10  (seawater).  The  eigen¬ 
values  are  from  the  database  described  above.  Similar  dependence  is  exhibited  by  the  modes  at 
other  frequencies  (e.g.,  10.2  and  11-1/3  kHz)  and  ground  conductivity  levels  (Ref.  22).  In  the 
nighttime  database  (see  Fig.  3.3-2);  the  mode  labeled  as  Mode  X  is  a  new  mode,  identified  in  re¬ 
cent  work  (Ref.  22).  Mode  X  is  found  to  exist  for  the  westerly  geomagnetic  bearings  and  low 
geomagnetic  latitudes  (typically  9m  <  15°)  at  all  Omega  frequencies  and  ground  conductivity 
levels.  Mode  X  eigenvalues  become  increasing  difficult  to  find  as  9m  increased  above  15°  and/ 
or  decreased  below  185°.  The  Mode  X  eigenvalues  are  shown  in  Fig.  3.3-l(b)  for  only 
9m  <,  10°.  Similar  Mode  X  behavior  is  seen  at  other  frequencies  and  ground  conductivity  lev¬ 
els.  The  causes  of  this  anomalous  Mode  X  behavior  are  currently  not  well  understood  and  need 
to  be  investigated.  In  the  nighttime  plot  (Fig.  3.3-1  (b)),  each  of  the  odd-numbered  modes  exhib¬ 
its  “trend  reversal”  in  its  9m  and  behavior.  Thus,  for  example  in  Fig.  3.3-1  (b),  over  the  en¬ 
tire  range  of  values,  the  /Jm  -dependence  of  Mode  1  for  9m  <,  5°  (and  Mode  3  for  9m  <  7°) 
is  opposite  in  shape  to  that  of  Mode  l’s  /?m  -dependence  for  9m  ^  6°  (and  Mode  3  for  9m  s  8°). 
Earlier  work  (Ref.  3)  did  not  identify  a  trend  reversal  in  the  nighttime  Mode  1  eigenvalues  with 
the  implication  that  the  /?m -dependence  for  the  eigenvalues  has  the  same  characteristic  shape 
for  9m  <  10°  as  that  found  for  9m  >  10°.  The  earlier  work  further  implies  that  higher-order 
odd  modes  have  a  similar  dependence  on  and  as  Mode  1.  The  nighttime  Mode  1  eigen¬ 
value  pattern  of  9m  and  dependence  implied  by  the  earlier  work  is  composed  of  portions  of 
Mode  1  and  Mode  X  as  defined  in  this  report.  For  example,  Mode  1  as  used  in  the  1980  PPC 
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Figure  3.3-l(a)  Eigenvalue  Plots  of  Daytime  Modes  1,2,3: 

Geomagnetic  Latitude/Geomagnetic  Bearing 


Geomagnetic  Bearing  (deg).  (3m 


model*  is  made  of  the  following  portions  of  the  eigenvalue  curves  in  Fig.  3.3-2:  (1)  Mode  1  for 
|<9m|  2:  6°,  (2)  Mode  1  for  |0m|  <  6°  and  90°  <  <  180°,  and  (3)  Mode  X  for  |0m|  <  5° 

and  180°  <  <  270°.  The  composite  behavior  of  this  mode,  referred  to  as  Mode  1:  Candi¬ 

date  n,  is  shown  in  Fig.  3.3-2(a).  The  other  possibility  is  that  the  Mode  1  eigenvalue  behavior  as 
a  function  of  and  6m  is  that  identified  in  Fig.  3.3-1  (b),  including  the  two  characteristic 
classes  of  shapes,  one  for  dm  >  6°  and  one  for  0m  <  6° .  The  composite  behavior  of  this  mode, 
referred  to  as  Mode  1:  Candidate  I,  is  depicted  in  Fig.  3.3-2(b). 

In  view  of  the  above  two  differing  candidate  structures  for  the  Mode  1  eigenvalues,  it  is 
recommended  that  the  model  structure  be  calibrated  using  functional  forms  derived  from  each 
of  the  two  candidate  Mode  1  structures.  If  sufficient  data  exists  with  the  appropriate  ranges  of 
0m  and  /?m ,  the  results  of  the  model  calibration  should  determine  which  candidate  best  explains 
the  observations.  Section  3.4  presents  the  functional  forms  for  the  two  candidates. 

Figure  3.3-3  depicts  the  0m  -  and  -dependence  of  the  phase  components  of  the  daytime 
Mode  1  for  conductivity  level  10  (seawater)  and  a  frequency  of  13.6  kHz.  Figures  3.3-4  and 
3.3-5  show  the  0m  and  /?m -dependence  of  the  computed  phase  components  of  the  nighttime 
Mode  1:  Candidate  I  and  Mode  1:  Candidate  II,  respectively,  for  the  same  ground  conductivity 
level  and  frequency  as  for  the  daytime  phase  components  plots.  The  detailed  expressions 
relating  the  signal  amplitude/phase  components  to  the  associated  mode  eigenvalue  and  signal 
frequency  are  given  in  Ref.  23.  Information  in  Figs.  3.3-3  through  3.3-5  is  derived  from  the  ei¬ 
genvalue  data  shown  in  Figs.  3.3-1  and  3.3-2.  The  phase  components  database,  as  described  in 
Section  3.4,  is  used  to  develop  the  structure  of  the  Mode  1  signal  phase  components  sub-model. 


3.4  SUB  MODEL  STRUCTURE 

This  section  presents  the  structure  of  the  geomagnetic/conductivity  sub-model  including 
two  candidates  (Mode  1:  Candidate  I  and  Mode  1:  Candidate  II)  for  the  nighttime  Mode  1.  To 
formulate  the  sub-model,  daytime/nighttime  plots  of  the  Mode  1  phase  components,  as  func¬ 
tions  of  frequency,  ground  conductivity  level,  geomagnetic  latitude,  and  geomagnetic  bearing 
were  examined  to  identify  signal  frequency  and  path  parameter-dependent  trends. 


*As  implied  by  the  phase  velocity  functional  forms  for  bearing/dip  angle. 


3-10 


Real  Fart  of  Eigenvalue  (deg) 


0-17317 

10-31-49 


Imaginary  Part  of  Eigenvalue  (deg) 


Figure  3.3-2(a)  Eigenvalue  Plots  of  Nighttime  Mode:  Candidate  II  as  Functions 
of  Geomagnetic  Latitude  and  Geomagnetic  Bearing 

0-17314 


Imaginary  Part  of  Eigenvalue  (deg) 

Figure  3.3*2(b)  Eigenvalue  Plots  of  Nighttime  Mode:  Candidate  l  as  Functions 
of  Geomagnetic  Latitude  and  Geomagnetic  Bearing 
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0  45  90  135  180  225  270  315  360 

Geomagnetic  Bearing  (deg) 

Figure  3.3-3(a)  Daytime  Mode  1  Phase  Velocity  Dependence  on  Geomagnetic 
Latitude/Geomagnetic  Bearing 
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G-173Z3 

10-31-89 


Figure  3.3-3(b)  Daytime  Mode  1  Excitation-Factor  Phase  Dependence  on  Geomagnetic 
Latitude/Geomagnetic  Bearing 


Examination  of  the  parameter  trends  in  the  plots  suggested  the  following  form  to  model 
the  phase  component  behavior  as  a  function  of  frequency,  conductivity-level,  and  illumination 
condition: 

Phase  component'!  jD/N 

for  a  specified  y  rt^r,J)/N  pD/N^  R  * 

combination  of  ^  H  rj 

fp,  O),  D/N  J  '**  (3.4-1) 

+  £  G^(@m,,3m)] 

i=  1 


where 

D/N  =  Day/Night 

fp  =  p*  frequency  (f,=10.2,  f2=ll  1/3,  and  f3=13.6  kHz) 

<7/  =  ground  conductivity  for  fh  conductivity  level  (/=  1,2 . 10) 

*  geomagnetic  latitude 
Pm  -  geomagnetic  bearing 
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Phase  Velocity  Variation  (cec/Mm) 
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10-31-89 


Figure  3.3-4(a)  Nighttime  Mode  1:  Candidate  I  Phase  Velocity  Dependence  on 
Geomagnetic  Latitude/Geomagnetic  Bearing 
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10-31-89 


Figure  3.3-4(b)  Nighttime  Mode  1:  Candidate  I  Excitation  Factor  Phase 

Dependence  on  Geomagnetic  Latitude/Geomagnetic  Bearing 
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Phase  Velocity  Variation(cec/Mm) 


G-17321 

10-31-89 


Geomagnetic  Bearing  (deg) 


Figure  3.3-5(a)  Nighttime  Mode  1:  Candidate  II  Phase  Velocity  Dependence 
on  Geomagnetic  Latitude/Geomagnetic  Bearing 


CJf’”a',D/N  =  jth  calibration  coefficient  associated  with  analytic  functional 

form  j  for  frequency  fp  and  conductivity  o/  for  day/night 

=  jth  analytic  functional  form  for  day  and  night  depending 
on  6m  and  pm 


=  i111  calibration  coefficient  associated  with  tabulated  functional 
form  i  for  frequency  fp  and  conductivity  o/  for  day/night 
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Excitation-Factor  Phase  (cec) 
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Geomagnetic  Bearing  (deg) 


Figure  3.3-5(b)  Nighttime  Mode  1:  Candidate  II  Excitation-Factor  Phase 

Dependence  on  Geomagnetic  Latitude/Geomagnetic  Bearing 
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Gfp’a'’N(0m,/?m)  =  i115  tabulated  nighttime  functional  form  for  frequency  fp  and 
conductivity  ct/  depending  on  Om  and 

J0^  =  Number  of  day  (D)/night  (N)  calibration  coefficients 

associated  with  the  analytic  functional  forms 

IN  =  Number  of  night  (N)  weighting  coefficients  associated  with 

the  tabulated  functional  forms. 

The  analytic  functional  form  is  a  combination  of  trigonometric  and/or  exponential  functions  of 
Bm  and  /?m,  whereas  the  tabulated  functional  form  is  a  look-up  table  of  data  values  for  each  con¬ 
ductivity  level  and  frequency  parametric  in  Bm  and  fim .  Note  that  the  analytic  functional  forms 
are  independent  of  frequency  and  conductivity  level,  whereas  the  tabulated  functional  forms  are 
specified  for  each  frequency  and  conductivity  level.  The  tabulated  functional  forms  describe  the 
nighttime  behavior  of  the  phase  components  for  each  of  the  two  nighttime  Mode  1  candidates 
(proposed  in  Section  3.3),  under  the  condition  6m  <  0m  ^  where  is  the  geomagnetic 

latitude  value  at  which  Mode  1  exhibits  reversal  of  its  characteristic  dependence  on  Bm  and  /Jm . 
Also,  varies  slightly  with  signal  frequency,  fp,  and  the  conductivity  level,  a/. 

The  sub-model  phase  components  dependence  on  Bm  and  should  be  determined  sep¬ 
arately  for  the  three  frequencies  and  five  of  the  ten  conductivity  levels.  The  reason  for  selecting 
five  of  the  ten  conductivity  levels  is  that  the  phase  component  value  for  any  of  the  missing  con¬ 
ductivity  levels  (i.e.,  levels  5,  6,  7,  8,  or  9;  see  Table  3.2-1)  can  be  accurately  obtained  by  linear¬ 
ly  interpolating  the  phase  component  values  between  the  conductivity  levels  4  and  10. 

The  methodology  for  developing  each  of  the  phase  components  of  the  sub-model  is  as 
follows.  First,  alternative  candidate  expressions,  consisting  of  different  sets  of  analytic/tabu¬ 
lated  functional  forms,  are  postulated  for  the  phase  component.  Second,  each  candidate  expres¬ 
sion  is  fit  to  the  data  (from  the  theoretical  database)  of  the  associated  phase  component,  using  a 
least-squares  estimation  technique.  Third,  the  rms  (root-mean-square)  fit  error  of  the  selected 
candidate  expressions  are  compared  along  with  the  number  of  calibration  coefficients  included 
in  the  corresponding  expressions.  Fourth,  and  finally,  based  on  a  tradeoff  between  minimizing 
the  number  of  calibration  coefficients  and  minimizing  the  fit  error,  a  “best-fit”  expression  is 
chosen  for  the  selected  phase  components.  Tables  3.4-1  and  3.4-2  present  the  calibration  coeffi¬ 
cients  and  the  analytic/tabulated  functional  forms  representing  the  selected  best-fit  expressions 
for  the  phase  velocity  and  excitation-factor-phase  components.  These  tables  show  the  theoreti¬ 
cal,  frequency-specific  values  of  the  calibration  coefficients  for  conductivity  level  10  (seawater). 
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Candidate  I  and  Mode  1:  Candidate  II 


Table  3.4-2  Excitation-Factor  Phase  Component*  Functional  Form; 
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late 


Note  that  the  calibration  coefficients  are  uniquely  defined  for  frequency  and  conductivity  level 
and  therefore  a  different  set  of  coefficients  will  be  obtained  for  different  combinations  of  fre¬ 
quency  and  conductivity  level.  Also,  the  tabulated  functional  forms  in  the  tables  will  be  different 
for  different  combinations  of  frequency  and  conductivity  level.  It  is  seen  that  the  analytic  func¬ 
tional  forms  are  same  for  both  nighttime  Mode  1  candidates  (Mode  1:  Candidate  I  and  Mode  2: 
Candidate  II);  however  the  tabulated  functional  forms  are  significantly  different  for  the  two 
candidates.  Although  Tables  3.4-1  and  3.4-2  include  a  single  tabulated  functional  form,  Gi,  for 
each  the  nighttime  phase  component  and  conductivity  level/frequency  combination,  it  would  be 
desirable  to  have  several  (instead  of  one)  tabulated  functional  forms  if  the  additional  calibration 
coefficients  associated  with  the  additional  tabulated  forms  can  be  adequately  calibrated  with  the 
available  observational  data. 

Table  3.4-3  gives  fit-error  performance  of  the  daytime  and  nighttime  phase  components 
for  conductivity  level  10  and  all  three  frequencies.  Since  the  nighttime  phase  components  for 
the  two  nighttime  Mode  1  candidates  differ  only  in  the  tabulated  functional  form  which  has 
no  fit  error,  the  performance  of  each  of  the  phase  components  is  the  same  for  the  two  Mode  1 
candidates. 


Table  3.4-3(a)  Performance  of  Best-Fit  Functional  Forms/Coefficients  for 
Phase  Velocity  Component:  Conductivity  Level  10 


ILLUMINATION 

NUMBER  OF 
CALIBRATION 
COEFFICIENTS 

RMS  (MAXIMUM)  FIT-ERROR  (cec/Mm) 
FOR  EACH  FREQUENCY  (kHz) 

10.2 

11-1/3 

13.6 

Day 

5 

.03  (0.07) 

.03  (.07) 

0.03  (0.05) 

Night 

5 

0.49  (1.48) 

0.57  (1.78) 

0.76  (2.41) 

Table  3.4-3(b)  Performance  of  Best-Fit  Functional  Forms/Coefficients  for 
Excitation-Factor  Phase  Component:  Conductivity  Level  10 


ILLUMINATION 

NUMBER  OF 
CALIBRATION 
COEFFICIENTS 

RMS  (MAXIMUM)  FIT-ERROR  (ccc) 

FOR  EACH  FREQUENCY  (kHz) 

10.2 

11-1/3 

13.6 

Day 

2 

0.05  (0.12) 

0.05  (0.14) 

0.05  (0.12) 

Night 

4 

0.21  (0.68) 

0.33  (0.89) 

0.76  (2.50) 
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3.5  SUMMARY 


The  geomagnetic/conductivity  sub-model  of  the  1990  PPC  model  characterizes  the  de¬ 
pendence  of  the  Omega  signal’s  Mode  1  phase  component  on  the  spatially-varying  geophysical 
parameters  (ground  conductivity  and  geomagnetic  field)  of  the  signal  path  for  day  and  night  so¬ 
lar  illumination  conditions,  and  signal  frequencies  of  10.2, 11-1/3, 13.6  kHz.  The  signal  phase  is 
described  in  terms  of  the  two  phase  components  (phase  velocity  and  excitation-factor  phase)  of 
the  Mode  1  signal.  Theoretical  phase  components  data  parametric  in  signal  frequency  and  spa¬ 
tial  path  parameters  for  both  day  and  night  conditions  are  obtained  using  a  waveguide-mode 
model  of  VLF  signal  propagation  prediction.  This  database  is  used  to  deduce  the  structure  of 
each  of  the  phase  components.  An  examination  of  the  plots  displaying  the  parametric  depen¬ 
dence  of  the  mode  eigenvalues  and  phase  components  indicates  that  the  nighttime  Mode  1  be¬ 
havior  (as  shown  in  Fig.  3.3-l(b))  is  significantly  different  from  that  implied  in  the  1980  PPC 
model.  Therefore,  the  sub-model  structure  is  developed  for  two  alternative  candidates  for  the 
nighttime  Mode  1,  referred  to  as  Mode  1:  Candidate  I  and  Mode  1:  Candidate  II.  Several  plausi¬ 
ble  candidate  structures  for  each  of  the  day  and  night  phase  components  are  postulated.  These 
candidates  are  calibrated  with  theoretical  data  and  a  best-fit  candidate  is  selected  for  the  sub¬ 
model.  It  is  recommended  that  this  structure,  as  a  part  of  the  1990  PPC  model,  be  calibrated 
with  observational  data  to:  (1)  select  the  correct  nighttime  mode  from  the  proposed  two  candi¬ 
dates,  and  (2)  determine  the  values  of  the  calibration  coefficients. 
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4.  OVERALL  PPC  MODEL  STRUCTURE 

4.1  INTRODUCTION 

The  temporal  and  spatial  sub-models,  introduced  in  Chapters  2  and  3,  describe  signal 
phase  variations  under  conditions  in  which  one  of  the  categories  of  parameters  (spatial  or  tem¬ 
poral)  varies  from  segment-to-segment  while  the  other  category  has  its  parameters  fixed  at 
nominal  values.  Thus,  in  the  diurnal  sub-model,  the  day-to-night  scaling  is  selected  as  that  ex¬ 
pected  for  a  mid-latitude,  north-south  path  over  seawater.  Similarly,  for  the  geomagnetic/lati¬ 
tude  sub-model,  the  local  ionosphere  is  considered  either  “day”  ( x  =  0)  or  night  (%  =  n).  The 
method  by  which  the  spatial  and  temporal  phase  variations  are  combined  to  provide  a  total 
phase  variation  is  specified  by  the  overall  model  structure.  This  chapter  explains  the  roles  of  the 
two  sub-models  in  calculating  the  total  phase  prediction/propagation  correction  (PPC).  Also,  a 
suggested  procedure  is  given  for  calibrating  the  model  coefficients  including  the  onset/recovery 
time  parameters.  Finally,  the  1980  (currently  used)  and  1990  (proposed)  PPC  model  structures 
are  compared. 

4.2  COMBINING  THE  SPATIAL/TEMPORAL  SUB-MODELS 

In  developing  the  diurnal  model,  the  signal  phase  variation  is  described  in  terms  of  the 
variation  in  effective  ionospheric  reflection  height  with  solar  zenith  angle.  For  the  geomagnetic/ 
conductivity  sub-model,  waveguide-mode  solutions  for  the  signal  phase  are  specified  for  a  ho¬ 
mogeneous  waveguide  with  an  upper  ionospheric  boundary,  at  a  fixed  reflection  height  (either 
day  or  night),  composed  of  a  plasma  magnetized  in  various  directions  with  respect  to  the  propa¬ 
gation  path  and  a  lower  boundary  comprising  a  dielectric  medium  with  various  electrical  con¬ 
ductivities.  In  terms  of  the  waveguide-mode  model,  the  phase  contributions  of  the  geomagnetic 
field  (upper  boundary)  and  phase  contributions  of  the  surface  conductivity  (lower  boundary) 
are  coupled,  i.e.,  they  cannot  be  independently  specified  in  computing  phase.  Moreover,  each 
component  functional  form  making  up  the  geomagnetic/conductivity  sub-model  is  uniquely 
scaled  to  represent  its  particular  contribution  to  the  phase  for  day  or  night.  Aside  from  the  scal¬ 
ing,  however,  the  evolution  of  the  component  spatial  functional  forms  between  the  day  and  night 
extremes  of  effective  reflection  height  should  be  the  same.  This  fundamental  assumption  ap¬ 
pears  to  be  consistent  with  both  theory  and  observation. 
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Thus,  each  spatial  functional  form  Fj  is  decoupled  from  the  evolutionary  (temporal) 
form  E(x)  and  the  resulting  local  space/time  phase  contribution  from  functional  form  i  is 

A^i  =  QFiEfr)  (4.2-1) 

where  Q  is  the  model  coefficient  for  spatial  functional  form  i  and  the  expression  for  E(x)  is  giv¬ 
en  by  either  Eq.  2.3-5  (daytime)  or  Eq.  2.4-6  (nighttime).  Essentially,  Eq.  4.2-1  states  that  func¬ 
tional  form  Fj  evolves  in  time  according  to  E(x)-  An  important  modification  to  Eq.  4.2-1  is 
made  depending  on  the  temporal  domain(s)  over  which  F  is  defined.  Equation  (4.2-1)  is  a  very 
general  expression  and  its  specific  form  depends  upon  the  temporal  domain  of  Fj  as  discussed 
below. 


Three  types  of  spatial  functional  forms  are  identified  which  evolve  uniquely  in  time  (i.e., 
solar  zenith  angle): 


1 .  Functional  forms  for  a  local  daytime  ionosphere  (0  <  x  <  ^thr) 

2.  Functional  forms  for  a  local  nighttime  ionosphere  (xthr  <  X  <  7t  ) 

3.  Functional  forms  defined  for  both  local  daytime  and  nighttime  iono¬ 
spheres  (0  <  x  <  »)• 


For  the  type  (1)  functional  forms,  phase  behavior  is  specified  for  daytime  (ionosphere 
under  direct  solar  control)  which  has  no  analog  at  night  (ionosphere  under  control  of  scattered 
Lyman-alpha  radiation).  Thus  as  the  solar  zenith  angle  increases  during  daytime,  this  type  of 
functional  form  must  contribute  less  to  the  total  predicted  phase  so  that  at  x  =  2Tthr  the  contri¬ 
bution  must  vanish.  This  means  that  the  daytime  evolution  function,  EoCf).  must  possess  the 
functional  dependence  on  %  given  by  Eq.  2.3-5  and  also  be  normalized  so  that  Ed(0)=1  and 
Ed(x)=0  for  %  ^  *THR-  Thus, 


c  /  \  _  Ioge(Ch(yrHR,  H^/ChQ;,  Hsc)) 
DW  logeChOtTK^-Hsc) 


(0  ^  X  *  xthr) 


(4.2-2) 


where  the  dependence  of  the  Chapman  Function  on  z'  has  been  suppressed.  Note  this  form  is  in¬ 
dependent  of  wavelength  since  the  spatial  functional  forms  (in  some  cases)  and  the  calibration 
coefficients  (in  all  cases)  will  be  wavelength-specific.  Figure  4.2-1  shows  a  schematic  plot  of 

EdCc). 


The  type  (2)  functional  forms  describe  phase  behavior  only  under  nighttime  conditions, 
i.e.,  phase  behavior  having  the  specific  functional  behavior  described  by  these  forms  is  not 
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Night-only  spatial  form:  A<j>  =  Cn  FnEn 


Spatial  forms  defined  for  both  night  and  day:  A<|>  =  (Cn-  Cd)  FEq  +  CdF 


Figure  4.2-1  Evolution  Functions  for  the  Three  Classes  of  Spatial  Functional  Forms 

theoretically  predicted  during  the  daytime.  Thus  as  the  solar  zenith  angle  decreases  during  the 
night,  a  smaller  contribution  is  expected  from  this  type  of  functional  form  until,  finally,  at 
X  -  Zthr.  the  contribution  vanishes.  The  nighttime  evolution  function  E^(x)  must  have  the 
functional  dependence  on  %  given  by  Eq.  2.4-6  and  also  be  normalized  so  that  En(jt)=1  and 
EnCc)=0  for  x  £  Xthr-  Hence, 

en(x)  =  ™  (ZThr  ^  X  &  *)  (4.2-3) 

1  +  cos  ZTHR 

As  for  the  daytime  case,  this  form  is  independent  of  wavelength  since  the  wavelength  depen¬ 
dence  is  absorbed  in  the  calibration  coefficients  Cn-  Figure  4.2-1  shows  a  schematic  plot  of 
EnOc)- 

Type  (3)  functional  forms  are  those  which  describe  specific  functional  behavior  of  the 
phase  predicted  for  both  daytime  and  nighttime  ionospheres.  The  night  and  day  model  coeffi¬ 
cients  associated  with  these  forms  furnish  a  measure  of  the  scale  of  the  day-to-night  phase  vari¬ 
ation  (diurnal  shift)  for  the  particular  functional  form.  Scaled  by  the  diurnal  shift,  the  phase 
variation  is  given  by  Eqs.  2.3-5  and  2.4-6  with  matching  at  %  =  ^thr.  A  default  night  variation 


4-3 


of  10%  of  the  total  diurnal  shift  is  based  on  preliminary  analysis  of  the  data  (Appendix  B).  Us¬ 
ing  these  conditions  and  normalizing  yields 


where 


EoCt) 


Hsc  loge  Ch(%,  Hsc) 
-  A  +  B 


X  <  Xthr 


Acosx  +  B 
-  A  +  B 


X  >  yfTHR 


Hsc  loge  Ch(ymR,  Hsc) 
(I-7XI  +  cos  ^thr) 


B 


Hsc  l°ge  Ch(*THR,  Hsc) 


/j  (1  + cos  Xthr)  j 
(1  -7)  (1  +  COS  XTHR ) 


(4.2-4) 


and  y  is  the  nighttime  phase  variation  expressed  as  a  fraction  of  the  diurnal  shift  (default 
y=0.1).  With  E0(x)  so  defined,  A<p  is  generalized  from  the  form  given  in  Eq.  4.2-1  to  explicitly 
show  how  the  difference  in  night  and  day  coefficients  determines  the  scale  of  the  diurnal  shift 
for  functional  form  Fj,  viz. 

A<Pi  =  (CN  -  CD)FiE0(*)  +  CDFi  (4.2-5) 

For  Fj  =  1 ,  the  purely  diurnal  phase  model  of  Chapter  2  is  recovered.  Figure  4.2-1  shows  a  plot 
of  E0(z). 


4.3  OVERALL  EXPRESSION  FOR  PREDICTED  PHASE  AND  PPC 

In  the  waveguide-mode  model,  the  predicted  VLF  signal  phase  measured  on  a  homoge¬ 
neous  path  from  a  transmitting  source  is 

<p  =  IcReV'  +  <Pe  +  0E  (4-3-1) 

where  k=l/vl  (A  =wavelength)  is  the  (assumed  constant)  wave  number,  Re  is  the  earth’s  radius, 
ip  is  the  radian  measure  of  the  (great-circle)  path,  and  <p£  are  the  excitation-factor  phases 
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at  the  transmitting  source  and  receive  point,  respectively.  For  an  inhomogeneous  path  (i.e.,  a 
real  path),  the  first  term  in  Eq.  4.3-1  is  replaced  by  a  path  integral  (or  sum).  For  the  path  seg¬ 
mentation  introduced  in  Chapter  2,  Eq.  4.3-1  becomes 

4>  =  XI  kijREAaj  +  +  *§  (4.3-2) 

*  i 

where  i  indexes  path  segment,  j  indexes  the  component  functional  form  of  the  spatial  sub¬ 
model,  and  Aaj  =  0.01  radian  for  all  i.  kjj  is  the  wave  number  corresponding  to  path  segment  i 
and  functional  form  j. 

In  order  to  express  the  phase  given  by  Eq.  4.3-1  in  terms  of  the  spatial  functional  forms/ 
model  coefficients,  a  relationship  is  needed  between  the  wave  number  and  the  model  coeffi¬ 
cient/functional  form  product.  Each  contribution  (j)  to  the  phase  velocity  (vj)  component  is 
actually  computed  as  (see,  e.g.,  Fig.  3.3-3(a),  vertical  scale) 

CjFj  =  (1-Vj/c)  j  -  (l-vj/cjko  (4.3-3) 

where  calibration  coefficient,  Cj,  has  units  of  cycles/Mm,  Fj  is  the  j,h  dimensionless  functional 
form  (and  has  magnitude  between  0  and  1),  ko  is  the  free  space  wave  number  in  units  of  Mm'1, 
and  c  is  the  speed  of  light  in  free  space.  The  wave  number  kj  is  related  to  the  phase  velocity  vj  as 


—  or  kjj  = 


_f_ 

Vij 


(4.3-4) 


where  f  is  the  frequency  and  the  i  is  attached  to  show  the  dependence  of  the  wave  number/phase 
velocity  on  path  segment.  Using  Eq.  4.3-3  to  obtain  Vj  (— ►  Vjj )  yields 


This  relation  together  with  Eq.  4.3-4,  produces  the  following  result  for  the  wave  number,  k,j 


kjj 


ko 

k, 


(4.3-5) 
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since  f  =  cko-  Inspection  of  Table  3.4-1  shows  that  |Cj|  <  0.20  cycle/Mm  and  |Fjj  |  <  1.  Thus, 
|  Cj  Fjj  |<  0.20  and  ko(10.2  kHz)  s  33  Mm'1  with  larger  values  for  the  other  frequencies.  Hence 
|Cj  Fjj/ko  |  «  1  and  only  first-order  terms  need  be  retained  in  the  binomial  series  expansion  of 
Eq.  4.3-5.  Thus 

CF 

kjj  =  ko(l  +  — J-1  +  ...)  s  ko+  CjFjj 
•to 

and  Eq.  4.3-2  becomes 


<P  =  koREa  +  ReXX  CjFjj  Aaj  +  +  <p%  (4.3-6) 

»  j 

If  Nj  is  defined  as  the  total  number  of  path  segments  and  Ns  is  the  number  of  path  seg¬ 
ments  to  the  point  of  first  ionospheric  interaction  (see  Fig.  2.5-4),  then  N-p-Ns  is  the  number  of 
segments  to  the  last  point  of  ionospheric  interaction.  Thus,  index  i  is  summed  from  Ns  to 
Nj-Ns  since  all  sub-model  components  are  defined  only  within  this  region  (recall  that  the  func¬ 
tional  dependence  of  phase  velocity  couples  ground  conductivity  and  ionospheric  parameters  in 
the  waveguide-mode  model).  The  so-called  “VLF  ground-wave”  effects  over  the  first  and  last 
Ns  segments  of  the  path  are  assumed  to  be  adequately  approximated  by  the  free-space  wave 
number.  The  excitation-factor  phase  components,  <f> <p£  are  assumed  to  depend  only  on  spa¬ 
tial  parameters  (magnetic  field,  path  orientation,  and  ground  conductivity)  at  the  segments  Ns 
and  Nj-Ns  since  these  are  the  first/last  point  of  ionospheric  interaction  and  thus  serve  to  define 
the  spatial  geophysical  environment  for  initial/final  excitation  of  the  wave  into/from  the  earth- 
ionosphere  waveguide.  Thus,  Eq.  4.3-6  may  be  written 

Nt- Nj  N,  Nj 

<p  =  koREa  +  Re  X  Z  cjpij  Aaj  +  X  X  ciFU  (4  3‘7) 

i  =  N$  j  =  1  i  =  NsNt-Ns  j  =  1 

where  the  unprimed  sums  are  associated  with  wave  number/phase  velocity  components  and  the 
primed  sums  with  excitation  factor  components  ( Nf,  Nf  is  the  number  of  functional  forms  for 
the  unprimed,  primed  component). 

The  expression  for  predicted  phase,  Eq.  4.3-7,  is  valid  only  at  an  instant  of  time.  Depen¬ 
dence  of  the  phase  prediction  on  solar  zenith  angle  is  implemented  by  inserting  the  three  types 
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of  evolution  functions  (day,  night,  and  day/night)  according  to  the  prescriptions,  Eq.  4.2-1  or 
Eq.  4.2-5.  Thus 

Nt-Ns 

<6  =  koREa  +  Re  2  2  CfF^cfe)  +  Xt^E,^ 

i  =  Ns  j  j 

,  (4.3-8) 

+  £  (Cf-C^FijEote)  +  £  CPF,  Aa,  +  2  I  cjF[j 

j  j  _  i  =  NsJ4i-Ns  j 

where  the  superscript  D(N)  on  the  functional  form  indicates  that  day-only  (night-only)  forms 
are  to  be  invoked;  no  superscript  implies  that  the  form  is  invoked  both  night  (with  coefficient 
CN)  and  day  (with  coefficient  CD).  The  primed  terms  again  refer  to  the  excitation  factor  phase 
component,  evaluated  at  segments  Ns  and  Nt-Ns,  which  evolves  with  solar  zenith  angle  in  the 
same  fashion  as  phase  velocity.  Since,  according  to  Table  3.4-2,  there  are  no  excitation  factor 
forms  defined  for  day  only*,  CD  FD'  =  0  and  the  excitation  factor  contribution  to  Eq.  4.3-8  be¬ 
comes 

I  XCjFii"  I  X  Cr'Fj'ENtoO  +  2<cf  -  Cf'^EoCCi)  +  2  Cf'Fi, 

i^NJVNs  j  i  =  Ns  J'J-j— Ns  j  j  j 

By  definition,  the  propagation  correction  (PPC)  is  (Ref.  3) 

PPC  ■  NOMINAL  PHASE  -  PREDICTED  PHASE 

where  the  nominal  phase  is  defined  for  a  great-circle  path  subtending  a  radians  as 

NOMINAL  PHASE  =  0.9974  koRE<z 

Thus,  the  PPC  becomes 

PPC  =  -0.0026  koREa-S 

where  S  is  that  part  of  the  predicted  phase  (Eq.  4.3-8)  which  involves  the  calibration  coeffi¬ 
cients,  i.e., 

S  =  0-  koREa 

•Note  that  Table  3.4-2  contains  day  and  night  functional  forms  for  the  excitation-factor 
phase  component.  However,  the  two  daytime  functional  forms  are  identical  to  two  of 
the  nighttime  forms  so  that  there  are  no  unique  daytime  forms. 
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4.4  MODEL  CALIBRATION  CONSIDERATIONS 


By  definition,  a  semi-empirical  model  contains  both  theoretical  and  observational  com¬ 
ponents.  Theory  is  used  to  determine  how  the  signal  phase  varies  with  a  single  parameter  (e.g., 
ground  conductivity)  holding  all  other  parameters  constant.  Analytical/tabular  approximations 
to  these  single-parameter  variations  comprise  the  functional  forms  described  in  Chapter  3.  To 
derive  a  complete  expression  for  phase  based  on  the  known  contributing  variables  is  intractable 
from  a  theoretical  viewpoint  (for  all  but  exactly  linear  models).  In  a  semi-empirical  formulation, 
calibration  coefficients  are  used  with  a  linearized  model  to  approximate  the  exact  dependence 
of  phase  on  the  dependent  variables.  Observational  data  is  then  enlisted  to  determine  the  best 
estimates  of  the  calibration  coefficients/parameters. 

The  calibration  procedure  must  be  carefully  chosen  due  to  the  parameter  coupling  in¬ 
herent  in  the  model  —  especially  the  spatial-temporal  coupling.  The  recommended  calibration 
procedure  addresses  the  three  different  categories  of  calibration  coefficients/parameters  ac¬ 
cording  to  the  following  sequence: 

1)  Daytime  coefficients  (both  day-only  and  day/night) 

2)  Nighttime  coefficients  (both  night-only  and  day/night) 

3)  Onset/recovery  time  parameters. 

This  sequence  (which  may  require  more  than  one  cycle)  specifies  calibration  of  the  daytime  co¬ 
efficients  first  because  these  coefficients  should  be  the  most  stable  and  well-defined.  Once 
these  are  determined  from  phase  data  over  the  appropriate  observation  hours  (as  defined  be¬ 
low),  the  nighttime  coefficients  are  computed  using  default  onset/recovery  time  parameter  val¬ 
ues.  Finally,  the  onset/recovery  time  parameters  are  determined  from  single-station  analog 
phase  data  (as  discussed  below). 

To  isolate  the  daytime  coefficients,  phase  data  is  selected  for  those  hours  producing  the 
minimum  predicted  phase  for  a  given  path.  Here,  “minimum”  is  used  loosely  since  the 
“midday”  phase  behavior  is  usually  quite  shallow  (see  Fig.  2.5-6).  To  define  this  shallowness 
more  quantitatively,  Eq.  4.2-2  may  be  used  to  show  that 

0.963  £  Ed(x)  £  1.000  for  0°  £  x  £  50° 

for  default  values  of  threshold  solar  zenith  angle  and  scale  height.  Moreover,  from  Eq.  4.2-4,  it 
can  be  shown  that 

0  £  E0(x)  £  0.0335  for  0°  £  %  £  50° 
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Thus  in  Eq.  4.3-8,  the  first  term  inside  the  brackets  (which  is  defined  for  night-only)  is  zero  and 
the  fourth  term  may  be  neglected  in  comparison  to  the  second  term  for  solar  zenith  angles  less 
than  50°*.  Thus,  Eq.  4.3-8  becomes 

Nr-Nsf  1 


0  =  koRga  +  Re  X 


I  CPfPEdCg)  +  2  CjFjj 


Act* 


(4.4-1) 


i  =  NsL  j 


j 


+  I  ZcP'Fi 

i  =  Ns^r-Ns  j 


0  <  x  <  50° 


Fora  given  path  and  time,  Eq.  4.4-1  can  be  used  with  observations  for  which  x  ^  50°  at  all  path 
segments  to  estimate  the  coefficients  CjD  and  Cj5’.  It  may  be  convenient  to  modify  this  procedure 
in  actual  practice  to  accept  phase  data  only  for  those  paths/hours  in  which  the  phase  change 
from  the  “minimum  day”  value  is  less  than  some  threshold  Y,alue  or  acceptable  fraction  of  the 
default  diurnal  shift. 

Once  the  daytime  calibration  coefficients  are  determined,  the  nighttime  coefficients  (CjN 
and  CjN)  are  estimated  with  the  use  of  path  phase  observations  measured  during  nighttime  con¬ 
ditions.  As  in  the  case  of  the  daytime  observations,  the  nighttime  measurement  hours  are  se¬ 
lected  so  that  the  appropriate  evolution  function  (En(^))  is  approximately  independent  of  x 
This  is  done  to  reduce  the  sensitivity  of  En  to  the  parameter  xth r  which  is  determined  in  the 
next  calibration  step.  Here,  an  iterative  approach  can  be  used  in  which  measurement  within  a 
given  range  of  x  are  used  to  calibrate  the  night  coefficients,  using  Eq.  4.3-8  and  a  default  value 
of  £thr-  ^  subsequent  calibration  yields  a  substantially  changed  value  of  ^thr.  then  the  night 
coefficients  should  be  recalibrated  using  measurements  from  a  smaller  range  of  x  •  Alternative¬ 
ly,  if  subsequent  calibration  of  xthr  yields  a  value  close  to  the  default  value,  then  the  night  coef¬ 
ficients  can  be  recalibrated  using  measurements  corresponding  to  a  larger  range  of  x 
Nighttime  evolution  functions  for  two  ranges  of  x>  computed  using  Eq.  4.2-3,  are  as  follows 
(*thr  =  96°) 

0.9327  <  En(z)  <  1  for  160°  <  *  <  180° 

0.6011  <  EnOc)  <  1  for  130°  <  x  <  180° 

The  same  considerations  apply  to  E0(x)  which  has  similar  numerical  values  for  the  two  ranges 
of  x  used  for  En(^)  above. 

“Small  excitation  effects  have  also  been  ignored. 
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For  nighttime  conditions,  Eq.  4.3-8  may  be  written 

Nr-Ns  f 

*  =  koREa  +  Re  I  2  C|*fNen  (*)  +  2(Cf*  -  CfjFijE,,  Of,) 

i  =  Ns  j  j 

+  ICfFli"Lai+  £  fzcH^'E^+^cN'.Cf'JFjjE^i)  (4.4-2) 

j  J  i  =  NsNT-N,L  j  j 

+  Zcf*i 

j 

Since  the  Cj3  and  CP  are  determined  in  the  first  calibration  step,  the  only  unknowns  (other  than 
the  onset/recovery  time  parameters  embedded  in  En,  E0)  on  the  right-hand  side  of  Eq.  4.4-2  are 
the  coefficients  C.N  and  CP'. 

The  three  principal  onset/recovery  time  parameters,  defined  in  Chapter  2,  are  given  as 
follows: 

1.  Distance  from  transmitter  (receiver)  to  point  of  first  (last)  ionospheric 
interaction  (AL) 

2.  Threshold  solar  zenith  angle  (^thr) 

3.  Ionospheric  response  time  (td>tn)- 

To  show  how  these  parameters  are  combined  to  give  the  sunrise/sunset  onset/recovery  time, 
consider  a  path  at  the  onset  of  sunrise  in  which  the  transmitter  is  east  of  the  receiver.  It  can  be 
easily  shown  that  other  situations  (i.e.,  sunrise  recovery,  sunset  onset/recovery,  and  paths  with 
different  east/west  orientations)  produce  the  same  basic  expression  for  the  appropriate  time 
condition.  The  onset  time  must  be  referenced  to  an  epoch  whose  occurrence  is  reasonably 
“close”  in  time.  A  convenient  choice  is  the  time  at  which  ground  sunrise  occurs  at  the  transmit¬ 
ter.  Thus,  the  relative  path  sunrise  onset  time,  for  this  case,  is  defined  as 

AT  =  (path  sunrise  onset  time  as  observed  on  phase  measurement) 

-  (time  at  which  x  =  tt/2  at  transmitter) 

This  onset  time  can  be  separated  into  three  components  as  follows: 

AT  *  ATi  +  AT2  +  AT3  =  T'(d  =  AL)-T(*  =  nf 2,  d  =  0) 
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where 

AT!  =  T(x  =  w/2,  d  =  AL)  -  T(*  =  */2,  d  =  0) 

AT2  =  T(*  =  rrHR,  d  =  AL)  -  T(x  =  7t/2,  d  =  AL) 

AT3  =  T'(d  =  AL)  -  T(x  =  3Tthr,  d  =  AL) 

In  the  above,  T(^,  d)  is  the  time  at  which  the  solar  zenith  angle  is  x  at  distance  d  along  the  path 
(measured  from  the  transmitter)  and  T'(d)  is  the  time  at  which  the  portion  of  the  ionosphere  at 
distance  d  on  the  path  becomes  ionized  as  a  result  of  incident  solar  radiation  at  solar  zenith 
angle  *thr-  Depending  on  the  values  of  j^thr  and  AL  and  the  orientation  of  the  path,  the  quan¬ 
tity  AT2  may  be  positive  or  negative  (see  Fig.  4.4-1).  From  the  above  definitions,  the  compo¬ 
nents  of  AT  may  be  written 

AT!  =  ^  ;  AT2  =  ^  ;  AT3  =  rD 

S  X 

where  S  is  the  rate  at  which  the  path  becomes  illuminated,  A^  =  xthr  -  7t/2,  x  *s  die  rate  of 
change  of  solar  zenith  angle,  and  tq  is  the  daytime  ionospheric  response  time  (see  Chapter  2). 


T(z-Jt/2,d«0)  T(x  =  XTHR-d!BAL)  T(z-7t/2,d  =  AL)  T’  (d  =  A  L ) 

Figure  4.4-1  Time  Sequence  Leading  to  Sunrise  Onset  from  Reference 
Time  Epoch  (Solar  Zenith  Angle,  x  *  *s  */2  at  Transmitter) 
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To  account  for  the  uncertainty  in  these  parameters  and  to  provide  for  comparison  with  observa¬ 
tions,  calibration  parameters  Cl,  Cz  ,  and  C*  are  inserted  as  follows: 

AT  =  CL  ^  +  Cx  &  +  Qr 
S  x 

Observations  of  path  sunrise  onset  time  can  thus  lead  to  estimates  of  the  calibration  pa¬ 
rameters  Cl,  Cx  ,  and  Q  since  the  quantities  S  and  x  near  the  time  of  onset  are  easily  com¬ 
puted.  Moreover,  with  the  use  of  default  values  of  AL,  xthr>  and  r ,  the  numerical  estimates  of 
the  calibration  coefficients  are  equivalent  to  estimates  of  the  true  values  of  AL ,  ^thr,  and  x . 
The  coefficients  Cl,  Cx  ,  CT  can  be  specified  separately  for  first  (near  transmitter)/last  (near 
receiver)  points  of  ionospheric  interaction  and  sunrise/sunset.  If  reciprocity  is,  as  expected,  ap¬ 
proximately  valid,  then  the  coefficients  should  be  the  same  for  the  first/last  point  of  ionospheric 
interaction  but  the  sunrise  and  sunset  values  are,  in  any  case,  likely  to  be  different. 

Once  the  onset/recovery  time  parameters  are  calibrated,  the  day  and  night  coefficient 
calibrations  can  be  refined  using  all  the  hourly  data  available.  This  will  permit  a  more  precise 
calibration  of  all  unknown  parameters  and  thereby  provide  more  accurate  predictions. 


4.5  COMPARISON  OF  1980  AND  1990  PPC  MODEL  STRUCTURE 

The  1980  and  1990  PPC  models  are  both  based  on  a  waveguide-mode  model  of  wave 
propagation  and  employ  a  path  segmentation  scheme  for  both  temporal  and  spatial  sub-model 
contributions  to  the  phase.  As  a  result,  the  models  are  similar  in  many  respects  but  there  are  im¬ 
portant  and  even  fundamental  differences  which  are  briefly  described  in  this  section. 

Table  4.5-1  summarizes  the  principal  distinctions  between  the  two  model  structures  by 
comparing  the  temporal  and  spatial  sub-model  features  and  the  overall  model  structure.  The 
following  three  sub-sections  amplify  the  comparison  between  the  two  model  structures  given  in 
the  figure. 

4.5.1  Diurnal  Sub-model  Comparison 

Both  models  rely  on  solar  zenith  angle  as  the  principal  time-dependent  variable  for  tem¬ 
poral  phase  variation.  The  1980  model  (Ref.  3)  uses  an  empirically-based  piecewise  linear 
(three  linear  forms  connecting  phase  variation  and  cos  x)  diurnal  function  which  prescribes  an 
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Table  4.5-1  Comparison  of  1980  (Current)  and  1990  (Proposed)  PPC  Model  Structure 


1980  MODEL  1990  MODEL 


Piecewise  linear  interpolation  function 

Day  model  based  on  Chapman 
theory 

Diurnal 

Transition  dynamics  based  on  linear 
electron  loss 

Night  model  based  on  Lyman- 
alpha  scattering 

Sunrise/sun  >et  dump  schedule 

Accurate  prediction  of  onset/re¬ 
covery  times  by  use  of  cali brat- 
able  parameters 

Correlated  analytic  magnetic  forms 

Coupled  magnetic/conductivity 
forms 

Geomagnetic/ 

Conductivity 

Uncoupled  conductivity  forms 

Mix  of  independent  analytic  & 
tabular  forms 

Possible  erroneous  Mode  1  behavior 
on  certain  azimuths  near  equator 

Mode  1/Mode  X  behavior  speci¬ 
fied  in  equatorial  region 

Interpolation  of  spatial  forms  between 
night/day  using  piecewise  linear  diur¬ 
nal  function 

Evolution  functions  describe  time 
behavior  of 
spatial  sub-models 

Overall 

Structure 

Calibration  using  all-day  or  all- 
night  paths 

Data  from  all  hours  used  for  cali¬ 
bration 

Day  data  constructed  using  previous 
calibration  parameters 

Cyclical  calibration  involving  day / 
night  spatial 

parameters  and  onset/recovery 
time  parameters 

unchanging  phase  variation  at  nig.it  (x  2  98°).  The  1990  model  employs  a  Chapman  theory- 
based  phase  variation  with  solar  zenith  angle  which  applies  during  the  period  when  solar  photo¬ 
ionization  controls  the  effective  reflection  height  of  the  ionosphere.  The  two  prediction  methods 
yield  similar  results  for  £  —  0  but  can  differ  markedly  for  larger  x 

For^  >70°  on  a  path  segment,  the  1980  model  invokes  a  “dynamic”  phase  variation  re¬ 
sulting  from  a  differential  equation  analogous  to  an  ionosphere  continuity  equation  with  a  linear 
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loss  term*.  Sunrise  and  sunset  schedules  of  electron  “dumping”  are  also  invoked  in  the  1980  to 
account  for  excess  ionization  by  minority  ionospheric  constituents.  The  temporal  functional 
form  used  in  the  1990  model  embodies  changes  from  a  simple  loge  sec  x  behavior  to  a  steeper, 
more  complex  form  (i.e.,  Chapman  Function  begins  accounting  for  a  curved  ionosphere)  for 
X  S  70°.  No  dump  schedules  are  included  in  the  1990  model  because  the  observed  phase 
variation  ascribed  to  anomalous  electron  dumps  does  not  warrant  the  complexity  of  their  inclu¬ 
sion  in  the  prediction  model.  The  1990  diurnal  sub-model  focuses  on  predicting  a  particular  set 
of  time  parameters  to  which  phase  prediction  error  is  especially  sensitive:  sunrise/sunset  onset/ 
recovery  times.  Parameters  are  identified  in  the  1990  sub-model  which  are  used  to  accurately 
compute  the  onset/recovery  times. 

As  mentioned  above,  the  1980  model  predicts  a  fixed  phase  variation  for^  £  98°,  i.e., 
no  change  in  predicted  nighttime  phase.  The  1990  model  includes  a  nighttime  variation  with  so¬ 
lar  zenith  angle  based  on  satellite  observations  of  scattered  Lyman-a  radiation.  The  1980  model 
also  contains  a  long-term  (year-dependent)  functional  form  which  is  being  investigated  under  a 
separate  effort  (Ref.  16).  The  work  described  in  this  report  does  not  address  a  long-term  phase 
variation. 

4.5.2  Spatial  Sub-model  Comparison 

The  1980  PPC  model  contains  separate,  uncoupled  spatial  sub-models  for:  (1)  geomag¬ 
netic  latitude/bearing,  (2)  ground  conductivity,  (3)  auroral  zone,  and  (4)  polar  regions.  The 
1990  PPC  model  contains  a  single  sub-model  for  geomagnetic  latitude/bearing  and  ground  con¬ 
ductivity  which  accounts  for  the  known  coupling  between  these  two  categories  of  spatial  param¬ 
eters.  Modification  of  the  auroral  and  polar  sub-models  is  not  addressed  in  the  work  described 
by  this  report  so  that  these  two  sub-models  are  currently  planned  for  inclusion  into  the  1990 
model  in  their  present  form.  The  1980  geomagnetic  latitude/bearing  sub-model  contains  func¬ 
tional  forms  which  are  highly  correlated  while  the  1990  geomagnetic/conductivity  sub-model 
contains  essentially  independent,  uncorrelated  functional  forms.  All  functional  forms  of  the 
1980  sub-models  are  analytic  forms  (including  the  constant,  unity)  whereas  the  1990  geomagne¬ 
tic/conductivity  sub-model  contains  a  mix  of  analytic  and  tabular  functional  forms  to  minimize 
the  complexity,  increase  the  accuracy,  and  facilitate  the  calibration  of  the  model.  The  1980  geo¬ 
magnetic/bearing  sub-model  specifies  functional  forms  on  certain  nighttime  low-latitude, 

’Current  models  of  D/E  region  ionization  dynamics  include  a  quadratic  loss  term 
as  indicated  in  Eq.  2.2-3. 


4-14 


westerly-directed  path  segments  which  may  implicitly  model  the  behavior  of  a  higher-order 
mode  (Mode  X  see  Chapter  3).  The  1990  geomagnetic/conductivity  sub-model  includes  two  sets 
of  candidate  functional  forms  which  permit  selection  of  the  correct  propagation  model(on  night¬ 
time  path  segments  within  the  appropriate  range  of  geomagnetic  bearing  and  latitude)  using  the 
observed  data. 


4.5.3  Overall  Model  Structure  Comparison 

Both  the  1980  and  1990  PPC  models  use  an  interpolative  scheme  to  combine  the  spatial 
and  temporal  sub-models  but  with  important  individual  differences.  The  1980  model  uses  an 
empirically-derived  diurnal  function  (modified  by  the  “dynamic”  diurnal  function  under  local 
transition  conditions)  which  is  piecewise-linear  in  cos  x  t0  interpolate  all  spatial  functional 
forms  between  their  daytime  and  nighttime  levels  (as  given  by  the  calibration  coefficients). 
Thus,  the  phase  variation  at  a  path  segment  (ignoring  transition  and  dump  schedule  contribu¬ 
tions)  is  given  by 

A0  -  ]£Cj)Fj+  ^  f(CjN  ~  Cj^Fj  (1980  model)  (4.5-1) 
j  j 


where  C^C*3)  are  the  calibration  coefficients  for  nighttime  (daytime)  corresponding  to  the  j1*1 
functional  form,  Fj ,  and  the  sum  is  over  all  spatial  functional  forms.  The  diurnal  function,  f,  is 
one  for  nighttime  conditions  and  decreases  through  two  types  of  linear  functions  of  cos  x  to  a 
value  of  zero  at  normal  illumination  (sun  directly  overhead).  Those  functional  forms  which  ap¬ 
ply  to  daytime  only  are  selected  by  taking  CjN  =  0.  Thus,  the  predicted  phase  variation  for 
daytime-only  forms  is  (from  Eq.  4.5-1) 

A<f>  =  X  -  f)Fj  (sum  over  day-only  forms) 


which  gradually  decreases  in  magnitude  as  x  increases  until  nighttime  (f=l)  is  achieved  at 
which  A<p  =  0 .  Similarly,  those  functional  forms  defined  for  night  only  are  selected  by  taking 
CjD  =  0.  In  this  case,  the  predicted  phase  variation  is  (from  Eq.  4.5-1) 


ZcffFj 


NIGHT  j 
ONLY  1 


(sum  over  night-only  forms) 


The  above  form  is  zero  only  for  f=0  (sun  directly  overhead)  but  for  all  other  values  ofx  it 
is  non-zero.  This  means  that  the  functional  forms  specified  for  nighttime-only  are  invoked  for 


virtually  all  local  daytime  conditions  (amplitude  of  the  contribution  increasing  with  increasing 
X)  as  well  as  nighttime  conditions. 


The  1990  PPC  model  avoids  this  difficulty  by  specifying  three  separate  evolution  func¬ 
tions  depending  on  the  temporal  domain  of  definition  for  the  functional  form:  (1)  daytime-only 
(2)  nighttime-only  and  (3)  daytime  and  nighttime.  In  each  case,  the  actual  form  of  the  variation 
is  based  on  either  the  Chapman  photoionization  model  (daytime)  or  scattered  Lyman-a  radi¬ 
ation  (nighttime)  but  the  forms  are  displaced/normalized  so  as  to  apply  only  in  the  appropriate 
temporal  domain. 

The  1980  PPC  model  is  structured  to  be  calibrated  using  data  recorded  only  on  signal 
paths  which  are  either  all-day  or  all-night.  Data  on  transition  paths  are  not  used  in  the  calibra¬ 
tion  since  no  calibration  coefficients  appear  in  the  diurnal  function  (or  “dynamic”  diurnal  func¬ 
tion).  This  means  that  large  observed  errors  on  transition  paths  will  not  be  corrected  by  the 
usual  process  of  coefficient  calibration.  Daytime  data  is  modified  (using  an  earlier  set  of  cali¬ 
bration  coefficients)  by  subtracting  the  time-dependent  component  of  the  predicted  phase  prior 
to  the  conventional  daytime  calibration  procedure.  Since  the  time-dependent  portion  depends 
on  the  old  calibration  coefficients,  a  cyclical  procedure  is  indicated  but  not  normally  carried  out. 
Because  the  diurnal  function,  f,  is  one  at  night,  the  nighttime  coefficient  calibration  requires  no 
a  priori  data  modification. 

The  1990  model  invokes  a  3-step  calibration  cycle  in  which  the  daytime  coefficients  are 
first  calibrated  from  data  on  paths  limited  to  a  certain  range  of  X-  In  the  second  step,  nighttime 
coefficients  are  also  calibrated  on  paths  limited  to  a  certain  range  of  x ■  The  limitation  on  x  in 
the  first  two  steps  is  primarily  due  to  the  need  to  reduce  the  sensitivity  of  the  terms  multiplying 
the  calibration  coefficients  to  the  onset/recovery  time  parameters  (in  particular,  ^thr)  as¬ 
sumed.  The  third  step  is  the  calibration  of  the  onset/recovery  time  parameters  using  analog  data 
from  which  sunrise/sunset  onset/recovery  times  are  extracted.  Comparison  of  the  calibrated  on¬ 
set/recovery  time  parameters  with  default  values  will  determine  the  revised  range  of  x  allowed 
for  paths  used  in  the  second  cycle  of  daytime/nighttime  coefficient  calibration.  Once  stable  val¬ 
ues  of  the  onset/recovery  time  parameters  are  obtained,  calibration  of  the  daytime  and  night¬ 
time  coefficients  can  be  performed  using  data  on  paths  with  any  range  of  solar  zenith  angle. 
This  means  that  a  well-balanced  set  of  data  on  all  paths  and  at  all  times  is  used  to  calibrate  the 
model. 
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4.6  SUMMARY 


The  model  structure  developed  in  this  report  specifies  a  means  for  combining  the  spatial 
sub-model  describing  the  phase  changes  due  to  local  variations  in  geomagnetic/conductivity  ef¬ 
fects  with  the  temporal  dependence  given  in  the  diurnal  sub-model.  The  functional  forms  for 
wave  number/phase  velocity  and  excitation-factor  phase  are  defined  for  two  temporal  extremes: 
day  (x  —  0)  and  night  (x  -  n ).  The  physics-based  diurnal  forms  (introduced  in  Chapter  2)  are 
shifted/normalized  (but  retaining  the  same  basic  form)  to  serve  as  evolution  functions  which 
modify  the  spatial  forms  at  intermediate  values  of  solar  zenith  angle.  Separate  evolution  func¬ 
tions  are  defined  for  three  types  of  spatial  functional  forms:  (1)  day-only,  (2)  night-only,  and  (3) 
day/night. 

Calibration  of  the  model  is  carried  out  in  a  three-step  cycle:  (1)  daytime  coefficient  cali¬ 
bration,  (2)  nighttime  coefficient  calibration,  and  (3)  onset/recovery  time  calibration.  Daytime 
coefficients  are  calibrated  using  data  recorded  at  hours  for  which  the  maximum  solar  zenith 
angle  is  less  than  about  50°  along  the  path.  These  coefficients  are  calibrated  first  because:  (1) 
the  predicted  phase  is,  for  the  indicated  range  of  solar  zenith  angles,  approximately  expressed 
in  terms  of  daytime  coefficients,  and  (2)  daytime  (midday)  phase  is  usually  inherently  stable 
and  repeatable.  Nighuime  coefficients  are  calibrated  next  using  data  recorded  at  hours  in  which 
the  path  solar  zenith  angles  are  in  some  appropriate  range  so  that  the  predicted  phase  is  essen¬ 
tially  independent  of  solar  zenith  angle  for  default  values  of  the  sunrise/sunset  onset/recovery 
time  parameters.  In  the  third  step,  the  onset/recovery  time  parameters,  which  include:  (1)  dis¬ 
tance  to  first/last  point  of  ionospheric  interaction,  (2)  threshold  solar  zenith  angle,  and  (3)  day / 
night  ionospheric  response  times,  are  calibrated  from  analog  phase  measurements  of 
onset/recovery  times.  Once  the  onset/recovery  time  parameters  are  established,  the  calibration 
cycle  can  be  repeated  without  limitation  on  the  path  range  of  solar  zenith  angle. 
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5. 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


The  development  and  structure  of  the  1990  PPC  model  are  summarized  in  the  first  sec¬ 
tion  of  this  chapter.  Conclusions  are  then  presented  regarding  the  physical  processes  addressed 
by  the  model,  spatial/temporal  coupling,  accommodation  of  calibration  data,  and  expected  ac¬ 
curacy.  The  1980  and  1990  PPC  model  structures  are  also  compared.  Finally,  recommendations 
are  made  regarding  sub-model  refinement,  additional  sub-model  development,  and  model  cali¬ 
bration  procedures/requirements. 

5.1  SUMMARY 

The  development  of  the  diurnal  sub-model  (presented  in  Chapter  2)  represents,  in  many 
ways,  a  new  approach  to  the  prediction  of  temporal  phase  by  focusing  on  physics-based  model¬ 
ing  of  the  ionosphere  and  on  those  predictive  features  to  which  the  phase  error  is  most  sensitive. 
The  importance  of  computing  correct  sunrise/sunset  onset/recovery  times  for  phase  prediction 
accuracy  is  demonstrated  to  motivate  the  analysis  of  parameters  which  determine  these  times. 
Background  information  on  ionospheric  physics  is  presented  to  provide  a  rramewcrk  for  the  de¬ 
velopment  of  daytime  and  nighttime  phase  variation  models.  The  daytime  phase  variation  is  de¬ 
veloped  using  a  Chapman  model  of  photoionization  dynamics  to  describe  how  the  effective  VLF 
wave  reflection  height  varies  with  solar  zenith  angle.  Effective  reflection  height  is  related  to 
phase  variation  to  obtain  the  functional  form  for  daytime  phase  variation.  The  nighttime  iono¬ 
sphere  is  not  controlled  directly  by  solar  radiation  so  that  other  sources  of  ionization  described 
in  the  literature  are  analyzed  in  terms  of  their  relative  dominance.  Solar  Lyman-alpha  radiation 
resonantly  scattered  from  the  earth’s  geocorona  is  the  dominant  nighttime  ionization  source  and 
satellite  measurements  of  Lyman-alpha  intensity  indicate  a  solar  zenith-angle  dependence  and 
sufficient  flux  to  ionize  the  nighttime  D/E-regions.  The  nighttime  phase  variation  is  determined 
by  using  an  approach  similar  to  Chapman’s  development  of  electron  production  in  the  iono¬ 
sphere.  Parameters  defining  sunrise/sunset  onset/recovery  times  are  described  in  terms  of  iono¬ 
spheric  properties  and  wave  propagation  features.  Results  of  a  path/time  phase  prediction 
computation  (using  only  diurnal  forms)  are  compared  with  observed  data. 

The  geomagnetic/conductivity  sub-model  is  a  fully  coupled  characterization  of  the  phase 
variation  due  to  spatial  variations  in  the  earth’s  magnetic  field,  signal  propagation  direction,  and 
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conductivity  of  the  earth’s  surface.  This  is  the  only  spatial  sub-model  treated  herein;  other  spa¬ 
tial  sub-models,  such  as  those  specifying  phase  variation  over  the  auroral  zones  and  polar  re¬ 
gions,  are  not  addressed  in  the  work  covered  by  this  report.  The  spatial  phase  variation  is 
separated  into  two  components:  (1)  a  path  length-dependent  part  known  as  the  wave  number  or 
phase  velocity  component  and  (2)  a  path  length-independent  part  referred  to  as  the  excitation 
factor  phase  component.  The  waveguide-mode  model  treats  the  total  signal  field  as  a  sum  of  in¬ 
dependent  signal  modes,  each  with  its  own  amplitude  and  phase.  By  design,  the  PPC  model  pre¬ 
dicts  only  the  phase  components  of  the  Mode  1  signal.  These  components  are  computed  from 
full-wave  waveguide-mode  models  in  terms  of  geomagnetic  latitude,  path  bearing  (with  respect 
to  geomagnetic  north),  and  ground  conductivity  for  both  day  and  night  for  each  of  the  three  fre¬ 
quencies,  10.2, 11-1/3,  and  13.6  kHz.  Functional  forms  are  developed  from  plots  of  the  theoreti¬ 
cal  phase  components  (as  a  function  of  each  spatial  parameter)  and  other  considerations.  For 
geomagnetic  latitudes  less  than  10  degrees,  phase  component  behavior  is  very  complex  so  that 
tabulated  functions  are  substituted  for  analytic  functional  forms.  Linear  combinations  of  func¬ 
tional  forms  are  fit  to  the  theoretical  phase  component  behavior  by  means  of  calibration  coeffi¬ 
cients.  At  low  geomagnetic  latitiudes  and  within  a  certain  range  of  geomagnetic  bearings,  the 
identity  of  the  Mode  1  signal  is  somewhat  obscure.  As  a  result,  two  candidate  forms  (tabular) 
are  developed  for  the  nighttime  phase  components  in  this  limited  region  of  parameter  space 
which  approximate  the  two  alternative  behaviors  for  Mode  1 .  Selection  of  the  correct  alternative 
is  possible  through  separate  calibration  of  each  candidate  functional  form  using  observational 
data. 


In  synthesizing  the  spatial  and  temporal  sub-models  into  an  overall  model  structure,  the 
basic  assumption  is  made  that,  although  the  diurnal  scaling  of  the  individual  spatial  functional 
forms  is  different,  their  evolution  in  changing  solar  zenith  angle  from  day  to  night  (or  night  to 
day)  has  the  same  functional  behavior  (given  by  the  evolution  functions)  for  all  spatial  function¬ 
al  forms.  Evolution  functions  are  developed  for  functional  forms  defined  within  three  temporal 
domains:  (1)  daytime-only,  (2)  nighttime-only,  and  (3)  both  day  and  night.  The  evolution  func¬ 
tion  for  a  given  temporal  domain  is  obtained  from  the  phase  behavior  developed  for  that  tempo¬ 
ral  domain.  The  expression  for  the  predicted  phase  and  PPC  is  written  using  the  defined 
evolution  functions  to  drive  the  solar  zenith  angle-dependence  of  the  functional  forms.  The 
model  structure  suggests  a  three-step  cyclical  calibration  procedure  in  which  daytime  coeffi¬ 
cients  are  first  calibrated  using  data  from  paths/times  having  an  allowed  range  of  solar  zenith 
angle.  In  the  second  step,  nighttime  coefficients  are  calibrated  using  data  from  paths/times 


5-2 


having  a  different  allowed  range  of  solar  zenith  angle.  The  third  step  involves  calibration  of  the 
onset/recovery  time  parameters  using  observed  onset/recovery  times  from  analog  data.  Once 
the  onset/recovery  times  are  established,  the  daytime  and  nighttime  coefficients  may  be  recali¬ 
brated  using  data  on  ail  path/times. 

Comparison  of  the  1980  and  1990  PPC  models  indicates  that,  for  the  temporal  depen¬ 
dence,  the  1990  model  daytime  phase  variation  follows  a  Chapman-type  behavior  with  solar  ze¬ 
nith  angle  up  to  the  day/night  threshold  value,  thus  excluding  any  special  contributions  from 
“transition”  effects  or  sunrise/suset  dump  schedules  included  in  the  1980  PPC  model.  The  1990 
model  predicted  nighttime  phase  varies  with  solar  zenith  angle  unlike  the  1980  model  in  which 
the  nighttime  phase  is  fixed.  The  phase  dependence  on  geomagnetic  field  (including  signal 
propagation  direction)  and  ground  conductivity  is  coupled  (as  predicted  by  theory) in  the  1990 
model  and  uncoupled  in  the  1980  model.  As  a  result  of  this  coupling,  the  spatial  functional 
forms  (all  analytic  for  the  1980  model  and  a  mix  of  analytic  and  tabular  for  the  1990  model)  for 
the  two  models  are  different.  Because  of  the  uncertain  identity  of  Mode  1  in  a  certain  range  of 
low  geomagnetic  latitudes  and  westerly  bearings,  the  1990  spatial  sub-model  contains  two  can¬ 
didate  forms  for  Mode  1  behavior  whereas  the  1980  model  assumes  that  only  one  of  the  forms 
characterizes  Mode  1  behavior.  The  1980  PPC  model  uses  a  diurnal  function,  piecewise  linear  in 
cosine  of  solar  zenith  angle,  as  a  means  of  combining  spatial  and  temporal  forms.  In  addition,  a 
“dynamic”  contribution  and  sunrise/sunset  “dumps”  are  added  during  transition  conditions. 
The  1990  model  uses  evolution  functions  for  spatial  functional  forms  defined  for  three  temporal 
domains  to  combine  spatial  and  temporal  forms.  Finally,  the  1980  model  is  structured  for  cali¬ 
bration  using  data  from  day-only  and  night-only  paths  whereas  the  1990  model  may  be  cali¬ 
brated  (via  a  three-step,  cyclical  procedure)  using  data  from  mixed  (day/night)  as  well  as 
day-only  and  night-only  paths. 

5.2  CONCLUSIONS 

The  structure  of  the  1990  PPC  model  presented  in  this  report  should  be  used  as  the  basis 
for  a  fully  calibrated  phase  prediction  model/algorithm  to  support  the  operational  needs  of 
Omega  users  worldwide.  The  features  of  this  semi-empirical  model  are: 

•  A  strong  adherence  to  physical  models  of  the  D/E-region  ionosphere  and 
VLF  signal  propagation 

•  Specific  focus  on  those  temporal/spatial  aspects  of  the  model  to  which 
phase  prediction  accuracy  is  most  sensitive 
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•  A  model  structure  which  is  amenable  to  calibration  using  all  path  data 
collected  on  a  24-hour  basis 

•  A  model  structure  which  guarantees  correct  temporal  invocation  of  ap¬ 
propriate  spatial  functional  forms 

•  A  minimal  number  of  independent  spatial  functional  forms  which  couple 
appropriate  geophysical  parameters  and  accurately  fit  theoretical  compu¬ 
tations. 

A  fully  operational,  properly  calibrated  model  based  on  the  structure  described  herein  is 
expected  to  provide  substantially  improved  phase  corrections.  This  expected  improvement  is 
based  on  a  comparison  of  the  1980  and  1990  PPC  models  in  terms  of  underlying  physical  mod¬ 
els,  expected  phase  prediction  accuracy  during  sunrise/sunset,  and  adaptibility  to  calibration  by 
available  observational  data.  The  potential  improvement  is  limited  by  the  fact  that  the  current 
auroral  and  spatial  sub-models  and  the  long-term  time  behavior  are  not  addressed  in  the  work 
covered  by  this  report.  Although  the  structure  of  auroral  and  polar  sub-models  is  not  analyzed, 
the  pertinent  geophysical  parameters  contributing  to  the  phase  variation  for  these  sub-models  is 
not  expected  to  couple  significantly  to  the  parameters  of  the  spatial  sub-model  considered  in 
this  report.  As  noted  earlier,  the  long-term  phase  dependence  is  not  addressed  in  the  work  cov¬ 
ered  by  this  report.  The  nature  of  the  long-term  phase  dependence  is  the  subject  of  a  separate  ef¬ 
fort  (Ref.  16)  and  its  effect  on  the  eventual  performance  of  the  model  is  not  known. 


5.3  RECOMMENDATIONS 

The  current  work  reported  herein  does  not  address  the  revision/modification  of  the  auro¬ 
ral/polar  sub-models.  It  is  recommended  that  these  sub-models  be  analyzed  in  light  of  current 
knowledge  of  polar/auroral  ionospheres  and  magnetosphere/ionosphere  interaction.  Upgrade  of 
these  semi-empirical  sub-models  to  a  synoptic-type  sub-model  (one  depending  on  near-real  time 
inputs)  should  also  be  investigated. 

The  1980  PPC  model  does  not  contain  an  equatorial  sub-model  (separate  contribution  to 
phase  variation  in  the  geomagnetic  equatorial  region).  However,  large  phase  prediction  errors 
observed  on  transequatorial  paths  suggest  the  need  for  such  a  sub-model.  It  is  recommended 
that  the  need  for,  and  the  physical  basis  of,  such  a  sub-model  be  investigated. 

The  1 1  1/3  kHz  calibration  coefficients  and  tabular  functional  forms  are  determined  only 
for  conductivity  level  10  since  a  separate  conductivity  sub-model,  per  se,  is  not  specifically 
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addressed  in  the  work  reported  herein.  It  is  recommended  that  the  1 1  1/3  kHz  calibration  coeffi¬ 
cients  and  tabular  functional  forms  for  the  geomagnetic/conductivity  sub-model  be  computed 
for  the  five  important  (in  terms  of  phase  variation  sensitivity)  conductivity  levels. 

Pending  a  decision  on  the  “other”  sub-models  (auroral,  polar,  equatorial,  and  long-term 
time  dependence),  it  is  recommended  that  the  1990  model  structure  presented  in  this  report  be  cali¬ 
brated  with  recent  years’  three-frequency  data  on  effectively  “one-way”  paths.  To  provide  observabil¬ 
ity  of  all  sub-model  functional  forms,  phase  difference  data  on  appropriate  paths  may  also  be 
required.  Following  calibration,  the  model  should  be  tested  and  compared  with  the  1980  PPC 
model  on  an  independent  database.  Finally,  the  model  calibration,  testing,  comparison  with  the 
1980  model,  and  operational  description  should  be  documented. 
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APPENDIX  A 

CHAPMAN  FUNCTION  APPROXIMATIONS 


The  Chapman  Function  is  important  in  the  quantitative  description  of  ionospheric  dy¬ 
namics.  Some  characteristics  and  an  approximate  analytic  expression  for  this  function  are  given 
in  this  appendix. 


The  Chapman  Function  evaluated  at  a  point  P  is  defined  (Ref.  7)  as 


Ch(x.  x) 
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where  r (h.'/)  is  the  optical  depth  at  point  P  which  is  at  altitude  h  above  the  surface  of  the  earth. 
X  is  the  corresponding  solar  zenith  angle  at  point  P,  a  is  the  absorption  cross  section  for  incident 
solar  radiation.  n{h)  is  the  neutral  density  at  altitude  h.  H  is  the  (assumed  constant)  scale  height. 
ds  is  a  differential  element  of  path  length  along  which  a  beam  of  radiation  travels,  and  .v  is  de¬ 
fined  by 

Re  t-  h 
x  =  - 

H 


where  Re  is  the  mean  radius  of  the  (assumed  spherical)  earth.  For  small  values  of  x  ,  the  curva¬ 
ture  of  the  ionosphere  (which  follows  the  curvature  of  the  earth)  can  be  safely  neglected  and  the 
Chapman  Function  is  adequately  approximated  by  sec  x  ■  However,  for  larger  values  of  solar  ze¬ 
nith  angle  satisfying 

x(l  —  sin^)  £  81  (A-2) 


the  sec  x  approximation  is  no  longer  valid  and  the  curvature  of  the  earth  must  be  taken  into  ac¬ 
count.  For  typical  values  encountered  in  VLF  signal  propagation  and  modeling  (Ref.  7),  x  s 
1290  and  the  solar  zenith  angle  corresponding  to  the  equality  in  Eq.  A-2  is  about  69.6°. 


Thus,  for  x  ~  69.6°,  Eq.  A-l  must  be  used  to  compute  the  Chapman  Function.  If  the 
neutral  density  is  assumed  to  decrease  exponentially  with  height  (exactly  true  for  an  ideal,  iso¬ 
thermal  gas  of  neutral  particles),  then  Eq.  A-l  can  be  approximately  integrated  (Ref.  7)  to  yield 


Ch(x,*) 


<  1  ±  erf 


X 


(A-3) 


A-l 


where  the  +  sign  is  used  for  x  >  tt/2  and  the  -  sign  is  used  for  x  <  jt/2. 


Equation  A- 3  is  the  fonn  used  for  computing  the  Chapman  Function  in  the  path/time 
computer  program  which  calculates  the  diurnal  phase  variation  for  VLF  signals  in  the  10-14 
kHz  range. 
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APPENDIX  B 

OBSERVATIONAL  PHASE  DATA 


This  appendix  describes  the  diurnal  characteristics  of  Omega  signal  phase  data  in  analog 
format  received  from  the  U.S.  Coast  Guard  Omega  Navigation  System  Center.  These  data  were 
analyzed  to  obtain  information  used  in  developing  and  testing  a  new  propagation  correction 
(PPC)  model. 

Section  B.2-1  presents  a  description  of  the  phase  data  analyzed  and  a  summary  of  the 
qualitative  diurnal  features  distilled  from  the  data.  Section  B.2-2  treats  a  selected  sample  of  ob¬ 
served  data  and  explains  certain  diurnal  features,  e.g.,  onset  and  recovery  times  of  phase 
changes  due  to  path  sunrise/sunset.  Sample  comparisons  of  observed  and  predicted  onset/re¬ 
covery  times  are  presented  in  Section  B.3.  This  section  also  describes  the  dependence  of  the 

T 

predicted  onset/recovery  times  on  specific  signal  propagation/ionospheric  parameters. 


B.l  PROCEDURES  AND  PRINCIPAL  FEATURES  OF  DATA  ANALYSIS 

Analog  phase  data,  recorded  on  effectively  “one-way”  Omega  signal  paths,  were  ex¬ 
amined  to  identify  the  fol’owing  time-dependent  features: 

•  Sunrise  onset/recovery  time 

•  Sunset  onset/recovery  time 

•  Qualitative  behavior  of  phase  on  all-day  paths 

•  Qualitative  behavior  of  phase  on  all-night  paths 

•  Time  of  minimum  phase  level  for  path  day 

•  Diurnal  shift 

•  Any  other  time-dependent  phase  behavior  characteristics. 

Some  of  this  data  is  used  to  develop  functional  forms  describing  time-dependent  phase  behavior 
for  the  1990  PPC  Diurnal  Sub-model.  Other  data  are  used  to  test/verify  phase  predictions  of  this 
sub-model. 


The  analog  data  analyzed  were  recorded  at  monitor  sites  co-located  with  Omega  trans¬ 
mitting  stations.  In  addition  to  the  phase  of  the  local  station  signal,  each  station  monitor  records 
the  phase  of: 

•  Six  other  station  signals  in  digital  format  (a  weak/poor  station  signal  is  cho¬ 
sen  as  a  calibration  channel) 

•  Four  other  station  signals  in  analog  format  (the  station  signals  used  in  sys¬ 
tem  synchronization) 

at  each  of  three  frequencies:  10.2,  11-1/3,  and  13.6  kHz.  Table  B.l-1  is  a  matrix  indicating  the 
paths  (defined  by  a  transmitting  station  and  receiver)  on  which  analog  data  is  recorded  at  the 
three  frequencies.  Table  B.l-2  gives  the  WGS  72/84  coordinates  (geodetic  latitude  and  longi¬ 
tude)  of  each  transmitting  station  and  its  associated  station  monitor.* 

Although  no  preciset  standard  time/frequency  signal  is  fed  directly  into  the  monitor  re¬ 
ceiver,  the  data  is  effectively  “one-way”  phase  since  the  phase  of  each  of  the  four  “remote” 


Table  B.l-1  Inter-station  Paths  (defined  by  transmitting  station  and 
receiver)  on  which  Analog  Data  is  Recorded 


TRANSMITTING 

STATION 

MONITOR 

RECEIVER 

A 

B 

c 

D 

E 

F 

G 

H 

A 

X 

X 

X 

X 

B 

X 

X 

X 

X 

C 

X 

X 

X 

X 

D 

X 

X 

X 

X 

E 

X 

X 

X 

X 

F 

X 

X 

X 

X 

G 

X 

X 

X 

X 

H 

X 

X 

X 

X 

*A  station  monitor  is  located  approximately  1-2  Omega  signal  wavelengths  from  the 
transmitting  antenna  which  is  sufficiently  far  so  that  the  radiated  signal  field  component 
dominates  the  electrostatic  and  inductive  field  components  and  close  enough  so  that  no 
diurnal  effect  is  seen  (no  appreciable  ionospheric  interaction). 

t“Precise”  means  an  accuracy  of  at  least  1  part  in  10n,  i.e.,  an  offset  of  about  one 
microsecond/day,  so  that  the  information  constrained  in  diurnal  traces  in  not  obscured 
by  oscillator  drift  problems. 
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station  signals  is  referenced  to  the  phase  of  the  local  station  signal,  which  is  essentially  constant. 
Thus,  the  phase  (modulo  one  cycle)  of  a  given  remote  station  signal  is  not  the  same  as  would  be 
obtained  by  a  receiver  time-synchronized  to  the  transmitted  signal;  the  difference  is  just  the  cu¬ 
mulative  phase  developed  between  the  local  transmitting  station  and  receiver  antennas.  Since 
this  difference  is  fixed  in  time,  however,  analysis  of  the  diurnal  behavior  (whrh  does  not  de¬ 
pend  on  absolute  phase  values)  is  unaffected.  The  effective  receiver  time  constant  is  adaptive, 
varying  between  1  and  5  minutes,  depending  on  the  tracking  error  in  the  phase-lock  loop. 

The  data  is  recorded  on  8-channel  chart  paper  with  pre-printed  hours  (15-20  minute  sub¬ 
divisions)  and  phase  values  between  0.00  and  1.00  cycle  (in  subdivisions  of  0.02  cycle  (2  cen- 
ticycles)).  Chan  speed  is  nominally  0.5  mm/minute  but  the  drive  motor  is  apparently  not 
precisely  controlled  so  that  adjustments  of  several  minutes  are  frequently  made  (and  annotated 
on  the  chan)  as  often  as  twice  per  day.  This  means  that  the  error  in  determining  specific  fea¬ 
tures,  such  as  onset  time,  is  limited  to  about  one  receiver  time  constant. 

The  analog  data  received  from  the  U.S.  Coast  Guard  Omega  Navigation  System  Center 
(ONSCEN)  consists  of  signal  phase  traces  on  the  paths  indicated  in  Table  B.  1-1,  at  the  three  fre¬ 
quencies  (10.2,  11-1/3,  13.6  kHz)for  certain  months  within  the  year  1988  (some  1987  and  1989 


Table  B.l-2  Transmitting  Stations  and  Station  Monitor  Coordinates 
for  WGS  1972/1984  Spheroids 


TRANSMITTING  STATIONS 

STATION  MONITORS 

ST  A  riON 

LATITUDE 

LONGITUDE 

LATITUDE 

LONGITUDE 

(DEGREES) 

(DEGP5ES) 

(DEGREES) 

(DEGREES) 

A 

66. 42017/66. 42019N 

013. 13681/013. 13697E 

66. 52930/66. 52932N 

012. 84530/012. 84545E 

B 

06. 30531/06. 30535N 

010. 66455/010. 66440W 

06. 42627/06. 42631N 

010.81398/C10.81382W 

r 

21. 40466/21. 40470N 

157. 83098/157. 83082W 

21. 52091/21. 52095N 

157. 99658/157. 99642W 

D 

46. 36591/46. 36594N 

098.33.  77/098. 33562W 

46. 55950/46. 55953N 

098. 63880/098. 63865W 

E 

20. 97418/20. 97414S 

055. 28974/055. 28990E 

20.91 194/20.9 1 19  IS 

055.5 1 695/055.5 1710E 

F 

43. 05358/43. 05355S 

065. 19093/065. 19078W 

42.75328V42.75325S 

065. 10094V065. 10079W 

G 

38. 48126/38. 48123S 

146. 93514/146. 93530E 

38. 39556/38. 39553S 

146. 66003/146. 66018E 

H 

34. 61470/34. 61474N 

129. 45349/129. 45365E 

34. 32471/34. 32474N 

_ i 

'Monitor  moved  in  November  1988;  previous  WGS  1972  coordinates  were: 
Latitude  =  43.22392°  South;  Longitude  =  65.27159°  West 
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data  were  also  included).  These  years  were  selected  to  obtain  representative,  recent  data  re¬ 
corded  during  periods  of  median  solar  activity.  Analysis  focused  on  10.2  and  13.6  kHz  data  in 
September,  since  this  month  includes  the  autumnal  equinox  which  is  expected  to  provide  a 
median  seasonal  component  of  phase  behavior.  The  days  selected  include  samples  near  the  be¬ 
ginning,  middle,  and  end  of  the  month.  A  reduced  sample  of  La  Reunion  data  was  analyzed 
since  this  station’s  annual  maintenance  is  scheduled  in  September. 


The  data  analysis  revealed  the  following  features  of  the  diurnal  phase  behavior: 


•  Due  to  external  noise  and  other  influences,  onset/recovery  times  have  a  max¬ 
imum  identification  accuracy  of  about  5  minutes;  paths  exhibiting  very  grad¬ 
ual  phase  onset/recovery  may  have  substantially  poorer  accuracy 

•  Path  sunset  onset  is  usually  a  more  gradual,  less  definite  process  than  path 
sunrise  onset 

•  The  value  of  the  phase  at  the  sunrise  recovery  time  generally  differs  from 
that  at  the  sunset  onset  time;  the  same  is  true  of  the  phase  at  sunset  recovery 
time  and  sunrise  onset  time 

•  Depressions  in  the  phase  following  sunrise  recovery  are  occasionally  found 
but  are  relatively  small  (less  than  10  cecs) 

•  Because  of  the  instabilities  observed  in  the  all-night  period  of  ihe  phase  trace 
on  certain  paths  (generally  attributed  to  modal  interference),  corresponding 
path  sunrise  onset  and  sunset  recovery  times  cannot  be  accurately  deter¬ 
mined  during  these  periods 

•  Sunrise/sunset  onset/recovery  times  are  difficult  to  identify  on  paths  exhibit¬ 
ing  a  double-diurnal  effect. 


The  double-diurnal  effect  referred  to  in  the  last  item  refers  to  a  characteristic  doubly  periodic 
phase  behavior  over  a  24-hour  period.  This  effect  occurs  in  phase  traces  on  paths  in  which  the 
path  vertex  (i.e.,  the  point  of  highest  latitude  on  the  path)  is  not  one  of  the  path  segment  end¬ 
points.  The  effect  is  generally  more  apparent  on  paths  in  which  the  latitude  of  the  path  vertex 
exceeds  the  latitude  of  the  terminator  vertex  (i.e.,  certain  polar  paths).  The  magnitude  of  the  ef¬ 
fect  is  diminished  by  the  fact  that  the  solar  zenith  angle  changes  relatively  little  over  24  hours  on 
polar  paths  since  the  path  is  “close”  to  the  earth’s  rotation  axis.  Figure  B.l-1  illustrates  the  dou¬ 
ble-diurnal  effect.  The  features  described  above  are  important  to  the  development  of  the  diurnal 
component  of  the  1990  t.'C  mode!  and  the  methods  used  in  its  calibration. 
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PATH  DEFINED  BY  TRANSMITTER  (T)  AND  RECEIVER  (R) 

Figure  B.l-1  Illustration  of  Double-Diumal  Effect;  North  Polar  View  during 
Northern  Hemisphere  Summer 


B.2  DESCRIPTION  OF  OBSERVED  DATA  SAMPLES 

Although  nearly  200  separate  phase  traces  were  analyzed  in  detail,  only  a  fraction  are 
presented  here  to  illustrate  the  behavior  of  real  phase  data  and  to  compare  predicted  and 
observed  sunrise/sunset  onset/recovery  times.  The  six  paths  selected  for  illustration  (three 
paths  plus  their  near-reciprocals)  have  large  north-south  components  to  minimize  the  spatial/ 
temporal  coupling  which  is  not  addressed  by  the  path/time  program  used  to  predict  the  onset/re¬ 
covery  times.  Signals  at  two  frequencies  (10.2  and  13.6  kHz)  are  displayed  for  each 
path/time-segment  to  permit  comparison. 

Figure  B.2-1  shows  the  observed  phase  traces  of  the  10.2  and  13.6  kHz  signals  on  the 
path  from  the  Liberia  transmitting  station  to  the  Norway  station  monitor.  Figure  B.2-1  (a)  illus¬ 
trates  the  first  12  hours  of  the  UT  day  /which  includes  sunrise  on  this  p?*' '  nd  Fig.  B.2-1  (b)  il¬ 
lustrates  the  next  12  hours  (including  path  sunset).  The  character..  .  .owl-shaped”  daytime 

‘The  paths  are  not  exactly  reciprocal  since  the  transmitting  antenna  and  monitor  are 
separated  by  20-50  km  at  a  station.  For  a  path  length  of  approximately  7,000  km, 
however,  the  difference  is  small,  so  that  the  paths  are  very  nearly  coincident. 
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G-17493 

11-8-89 


a)  0000-1200  UT 


b)  1200-2400  UT 


Figure  B.2-1  10.2  and  13.6  kHz  Phase  Data  from  the  Liberia  Transmitter  to 
the  Norway  Station  Monitor  on  26  September  1988 


trace  is  noted  between  the  hours  of  about  0600  and  1800  UT  and  a  slight  peak  is  observed  on  the 
nighttime  trace  near  2400  UT.  A  slight  post-sunrise  depression  occurs  on  both  10.2  and  13.6 
kHz.  The  arrows  in  the  figures  indicate  the  predicted  onset/recovery  times.  For  example,  in 
Fig.  B.2-1  (a)  the  arrow  prior  to  0600  UT  designates  sunrise  onset  and  the  arrow  following  0600 
UT  denotes  sunrise  recovery.  Similarly,  the  arrows  in  Fig.  B.2-1  (b)  indicate  sunset  onset  and  re¬ 
covery  times.  Note  that  sunset  onset  time  is  not  so  readily  apparent  as  the  other  times;  closer  ex¬ 
amination  shows  that  this  time  separates  the  bowl-shaped  daytime  interval  from  the  more  linear 
path  sunset  period. 

Figure  B.2-2  displays  the  path  reciprocal*  to  that  described  above,  i.e.,  from  the  Norway 
transmitting  station  to  the  station  monitor  at  Liberia.  The  phase  behavior  is  similar  to  that  of 
Fig.  B.2-1  except  for  greater  small-scale  instability  due  probably  to  the  greater  local  noise  at  the 
Liberia  monitor.  Note  that  the  13.6  kHz  trace  shows  somewhat  less  irregularity  than  that  for 
10.2  kHz,  likely  because  of  the  expected  higher  signal  strength  at  13.6  kHz. 

The  signal  phase  on  the  path  from  the  Argentina  transmitting  station  to  the  station  moni¬ 
tor  at  North  Dakota  is  shown  in  Fig.  B.2-3.  The  phase  trace  during  the  first  12  hours  of  the  UT 
day  shown  in  Fig.  B.2-3  (a)  is  highly  oscillatory  —  the  result  of  nighttime  modal  interference, 
characteristic  of  nighttime  transequatorial  paths  having  a  westerly  propagating  component. 
Note  the  relatively  fewer  oscillations  on  the  13.6  kHz  trace  than  the  10.2  kHz  trace  prior  to  0600 
UT,  corresponding  to  less  modal  interference  at  this  frequency.  Onset/recovery  times  are  shown 
in  the  figure,  but  are  not  particularly  meaningful  during  periods  in  which  modal  interference  is 
this  severe.  Note  especially  the  single  large  oscillation  following  the  predicted  sunrise  onset 
time.  This  oscillation  probably  results  from  modal  conversion  at  the  terminator  and  is  not  ob¬ 
served  at  sunset  recovery  —  possibly  because  of  the  different  path-terminator  orientations  at 
sunrise  and  sunset.  Sunrise  recovery  at  both  frequencies,  however,  is  cleanly  predicted  Fig¬ 
ure  B.2-3(b)  shows  primarily  the  daytime  phase  behavior  which  is  much  more  stable  than  the 
nighttime  phase  behavior.  This  figure  shows  the  asymmetry  between  sunrise  recovery  (begin¬ 
ning  of  all-day  path)  and  sunset  onset  (end  of  all-day  path). 

Figure  B.2-4  displays  the  reciprocal  path,  i.e.,  from  the  North  Dakota  transmitting  anten¬ 
na  to  the  station  monitor  at  Argentina.  It  is  immediately  evident  from  Fig.  B.2-4(a)  that  little  or 
no  modal  interference  is  observed  on  the  phase  of  the  signal  propagated  on  the  same  path  in  the 
opposite  direction  (geomagnetic  south  and  east).  Nighttime  and  daytime  (see  Fig.  B.2-4(b))  be- 
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G-17494 

11-8-39 


a)  0000*1200  UT 


b)  1200-2400  UT 


Figure  B.2-2  10.2  and  13.6  kHz  Phase  Data  from  the  Norway  Transmitter  to 
the  Liberia  Station  Monitor  on  25  September  1988 


(cecs)  _  (cecs)  _  (cecs)  _  (cecs) 


havior  is  quite  regular  and  the  predicted  onset/recovery  times  are  consistent  with  the  observa¬ 
tions. 


Phase  behavior  on  the  path  from  the  Japan  transmitting  antenna  to  the  station  monitor  at 
Australia  is  shown  in  Fig.  B.2-5.  The  daytime  and  sunset  periods  of  the  phase  traces  are  given  in 
Fig.  B.2-5  (a).  A  distinguishing  feature  of  the  trace  is  the  Sudden  Ionospheric  Disturbance  (SID) 
occurring  shortly  after  0000  UT  on  both  frequencies.  The  daytime  and  sunset  phase  behavior  is 
otherwise  quite  normal.  Fig.  B.2-5 (b)  displays  the  nighttime  and  sunrise  periods  of  phase  be¬ 
havior.  The  nighttime  phase  has  a  broad  maximum  between  1200  and  1800  UT. 

Figure  B.2-6  illustrates  the  phase  behavior  of  the  path  from  the  Australia  transmitting 
antenna  to  the  station  monitor  in  Japan,  which  is  the  reciprocal  of  the  previously-described  path. 
The  SID  noted  earlier  is  also  seen  in  Fig.  B.2-6(a).  The  propagated  signal  on  this  transequatorial 
path  has  a  westerly  geomagnetic  component  and  thus  would  be  expected  to  exhibit  modal  inter¬ 
ference.  Figure  B.2-6(b)  shows  a  distinctive  modal  signature  on  the  10.2  kHz  trace,  especially 
near  1800  UT.  Following  1800  UT,  the  phase  trace  exhibits  a  near  cycle  slip,  recovering  (in  the 
correct  cycle)  just  prior  to  sunrise  recovery.  In  contrast,  the  13.6  kHz  phase  trace  is  reasonably 
well-behaved,  although  a  slight  step  is  seen  during  path  sunrise. 


B.3  COMPARISON  OF  OBSERVED  AND  PREDICTED  ONSET/RECOVERY  TIMES 

Table  B.3-1  provides  a  comparison  of  observed  and  predicted  sunrise/sunset  onset/re¬ 
covery  times  for  the  phase  data  described  in  Section  B.2.  The  table  provides  a  sample  compari¬ 
son  only  and  should  not  be  construed  as  an  indicator  of  ultimate  prediction  accuracy. 

Onset/recovery  time  prediction  accuracy  is  difficult  to  measure  because  of  the  uncertain¬ 
ty  inherent  in  identifying  time  of  onset/recovery  using  the  analog  phase  trace.  In  some  cases,  the 
identification  is  obscured  by  modal  interference  or  noise.  In  other  cases,  however,  identification 
problems  arise  because  of  no  apparent  change  in  phase  behavior  near  the  expected  time. 
Identification  errors  also  result  from  the  assumed  symmetry  in  the  phase  at  sunrise  recovery 
(beginning  of  all-day  path)  and  sunset  onset  (end  of  all-day  path).  Both  of  these  problems  con¬ 
tribute  to  the  large  differences  between  observed  and  predicted  sunset  onset  times  shown  in 
Table  B.3-1.  Comparison  of  predictions  and  observations,  in  addition  to  analysis  of  phase  data 
on  diverse  types  of  paths,  greatly  improves  identification  accuracy. 

Prediction  of  onset/recovery  time  is  based  on  three  principal  elements: 
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a)  0000-1200  UT 
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13.6  kHz 


b)  1200-2400  UT 


Figure  B.2-5  10.2  and  13.6  kHz  Phase  Data  from  the  Japan  Transmitter  to 
the  Australia  Station  Monitor  on  25  September  1988 
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G-17498 

11-8-89 


13.6  kHz 


10.2  kHz 


a)  0000-1200  UT 


13.6  kHz 


b)  1200-2400  UT 


10.2  kHz 


Figure  B.2-6  10.2  and  13.6  kHz  Phase  Data  from  the  Australia  Transmitter  to 
the  Japan  Station  Monitor  on  25  September  1988 


Table  B.3-1  Observed  and  Predicted  Onset/Recovery  Times  for  Sample  Data 


FRQ. 

(kHz) 

SUNRISE  TIMES  (UT) 
OBSERVED/PREDICTED 

SUNSET  TIMES  (UT) 
OBSERVED/PREDICTED 

PATH* 

DATE 

ONSET 

RECOVERY 

ONSET 

RECOVERY 

B-A 

191 

26  SEP  88 

0500/0454 

0620/0612 

<N 
*— < 

oo 

T—i 

• 

• 

1900/1900 

B-A 

w 

26  SEP  88 

0455/0454 

0620/0612 

••  / 1 8 12 

1904/1900 

A-B 

mSM 

25  SEP  88 

0454/0448 

0615/0612 

1735/1812 

1858/1900 

A-B 

13.6 

25  SEP  88 

0448/0448 

0604/0612 

1800/1812 

1857/1900 

F-D 

10.2 

26  SEP  88 

**  /0954 

1156/1154 

2305/2306 

**  /0048 

F-D 

13.6 

26  SEP  88 

**  /0954 

1155/1154 

2225/2306 

0200/0048 

D-F 

10.2 

26  SEP  88 

0950/1000 

1200/1154 

2315/2312 

0102/0048 

D-F 

13.6 

26  SEP  88 

0940/1000 

1157/1154 

2312/2312 

0043/0048 

H-G 

10.2 

25  SEP  88 

1930/1948 

2100/2048 

0825/0854 

1005/0942 

H-G 

13.6 

25  SEP  88 

1938/1948 

2104/2048 

0815/0854 

0937/0942 

G-H 

10.2 

25  SEP  88 

**  /1948 

2102/2048 

0755/0854 

0920/0942 

G-H 

13.6 

25  SEP  88 

_ 

1948/1948 

2100/2048 

0800/0854 

1010/0942 

•  First  letter  indicates  transmitter;  second  letter  indicates  station  monitor  receiver 


*  *  Measurement  not  recorded  due  to  large  uncertainty  in  identifying  time 


•  Point  of  first/last  interaction  between  the  signal  wave  field  and  the  iono¬ 
sphere 

•  Solar  zenith  angle  threshold  separating  direct  solar  control  of  the  local  iono¬ 
sphere  from  control  by  nighttime  sources 

•  Ionization  response  time  of  the  local  ionosphere  following  exposure  to  an 
ionization  source. 

These  parame^rs  can  be  varied  to  match  the  observed  on^et/recovery  times.  The  predic¬ 
tions  given  in  Table  B.3-1  are  based  on  default  values  of  these  parameters  and  thus  are  not  “op¬ 
timized.”  Also,  the  predictions  are  made  at  6-minute  intervals  which  are  about  1-3  receiver  time 
constants  and  comparable  to  the  uncertainty  in  the  time  information  derived  from  the  analog 
charts. 


Figures  B.3-1  through  B.3-3  show  the  sensitivity  of  onset/recovery  times  to  each  of  the 
three  onset/recovery  time  parameters  listed  above.  The  effect  of  varying  the  first/last  point  of 
ionospheric  interaction  is  illustrated  in  Fig.  B.3-1.  Two  traces  are  shown  which  differ  in  effec¬ 
tive  ionospheric  path  length  by  2  x  0.14  radian  =  0.28  radian  (about  20%  of  the  total  path 
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Tkna  (Hours) 


Lagond: 

-  First/Last  point  of  interaction  »  0.15  radian  from  path  arid-points 

- First/Last  point  at  interaction  -  0.01  radian  from  path  end-points 

SSO  Sunset  Onaot 

SSR  Sunset  Recovery 


Figure  B.3-1  Illustration  of  Effect  of  First/Last  Point  of  Wave/Ionosphere 
Interaction  on  Diurnal  Phase  Behavior 

length).  Increasing  the  effective  ionospheric  path  length  on  a  path  of  fixed  total  length  leads  to 
earlier  onset  times  and  later  recovery  times  as  shown  in  the  figure.  Figure  B.3-2  illustrates  the 
effect  of  varying  solar  zenith  angle  threshold  on  diurnal  phase  behavior.  The  two  traces  dis¬ 
played  represent  probable  extremes  in  physically  realizable  solar  zenith  angle  thresholds  (91° 
and  100°).  Larger  solar  zenith  angle  thresholds  tend  to  shift  both  onset  and  recovery  to  later 
times  as  compared  to  smaller  thresholds.  In  this  case  the  effective  ionospheric  path  length  does 
not  change  in  varying  the  thresholds,  so  that  the  diurnal  shift  is  the  same.  Note  the  difference  in 
phase  level  of  the  two  traces  in  the  all-day  portion  of  the  path;  this  is  due  to  the  dependence  of 
the  daytime  evolution  function  on  the  solar  zenith  angle  threshold.  The  effect  of  the  different 
daytime  and  nighttime  ionospheric  response  times  is  shown  in  Fig.  B.3-3.  The  solid  trace 
corresponds  to  an  equal  daytime  and  nighttime  response  time  of  six  minutes  while  the  dashed 
curve  traces  the  diurnal  behavior  for  the  same  daytime  response  time  but  a  nighttime  response 
of  66  minutes  (Ref.  16).  This  day/night  response  time  difference  has  the  effect  of  delaying  the 
sunset  onset  and  recovery  times  (as  shown)  but  little  effect  on  the  corresponding  times  for  sun¬ 
rise  since  the  strong  solar  photoionization  takes  over  rapidly  from  the  weaker  nighttime  sources. 
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Total  Phase  -  297  Cycles  (CenUcyclei)  Total  Phase  -  247  Cycto*  (Centtcyctes) 


G-17J0! 
1 1-1-19 


Tim*  (Hour*) 

legend: 

-  Solar  Zenith  Anglo  Threshold  -  100* 

- Solar  Zenith  Angle  Threshold  »  91  * 

SSO  Sunset  Onset 
SSR  Sunset  Recovery 

Figure  B.3-2  Illustration  of  the  Effect  of  Solar  Zenith  Angle  Threshold  on 
Diurnal  Phase  Behavior 

G-17500 

U-«-«9 


OMSTA  Australia  at  HAWAI  Monitor  (10.2  kHz)  05  Sop  88 


Time  (Hours) 

legend: 

—  Dayflme/Nightime  Response  ■  8  Minutes 

- Daytime  Response  -  6  Minutes  ;  Nightime  Response  »  66  Minutes 

SSO  Sunset  Onset 
SSR  Sunset  Recovery 

Figure  B.3-3  Illustration  of  the  Effect  of  Different  Daytime  and  Nighttime 
Ionospheric  Response  Times  on  Diurnal  Phase  Behavior 
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